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Measurements with a one-electron quantum cyclotron determine the electron magnetic mo- 
ment, given by g/2 = 1.001 159 652180 73(28) [0.28 ppt], and the fine structure constant, a~ — 
137.035 999 084 (51) [0.37 ppb]. Brief announcements of these measurements [HE] are supplemented 
here with a more complete description of the one-electron quantum cyclotron and the new measure- 
ment methods, a discussion of the cavity control of the radiation field, a summary of the analysis 
of the measurements, and a fuller discussion of the uncertainties. 



I. INTRODUCTION 

A. The Electron Magnetic Moment 

Measurements of the electron magnetic moment \i 
probe the interaction of the electron with the fluctuat- 
ing vacuum, allow the highest accuracy determination of 
the fine structure constant, and sensitively test quantum 
electrodynamics (QED). For an eigenstate of spin S, 
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where g/2 is the magnitude of fi scaled by the Bohr mag- 
neton \ib — eh/(2m). 

For angular momentum arising from orbital motion, 
g/2 depends on the relative distribution of charge and 
mass and equals 1/2 if they coincide, for example cy- 
clotron motion in a magnetic field. For a point particle 
in a renormalizable Dirac description, g/2 = 1, and de- 
viations from this value probe a particle's interactions 
with the vacuum as well as the nature of the particle it- 
self, as with the proton whose g/2 « 2.8 arises from its 
quark-gluon composition. 



B. New Measurements of the Electron Moment 

Our new measurements, announced in 2006 [T] and 
2008 0, used a one-electron quantum cyclotron [3] to 
determine the electron g/2 to a 0.76 ppt and then to a 
0.28 ppt accuracy. The latter result, 

ff /2 = 1.001 159 652180 73(28) [0.28 ppt], (2) 

has an uncertainty that is 2.7 and 15 times smaller than 
the 2006 and 1987 measurements (Fig.[T]), the latter being 
a measurement that stood for nearly twenty years [4]. 
The electron g is measured with an uncertainty that is 
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FIG. 1. Measurements 1, 2,4 of the dimensionless magnetic 
moment of the electron, g/2, which is the electron magnetic 
moment in Bohr magnetons. 



2300 times smaller than has been achieved for the heavier 
muon lepton [5]. 

The central feature of the new measurements is the 
quantum jump spectroscopy of completely resolved cy- 
clotron and spin levels of a one-electron quantum cy- 
clotron [3J. A number of new methods were introduced 
to make this possible. 

1. A cylindrical Penning trap cavity that was in- 
vented for these experiments [6] imposes boundary 
conditions upon the radiation field as well as 
providing an electrostatic quadrupole potential 
in which a single particle can be suspended and 
observed [7]. 

2. The resulting cavity-inhibited spontaneous emis- 
sion, at a rate 10 to 50 times below the radiation 
rate in free space, gives the averaging time re- 
quired to resolve one-quantum transitions that 
are made when all detection systems are turned off. 

3. Stored electron plasmas [8HI0] and the damping of 
a single electron in this cavity [5] are used together 
to determine cavity frequency shifts and eliminate 
cavity shifts as a major uncertainty [2]. 

4. Blackbody photons that would cause unwanted 
quantum jumps are eliminated by lowering the 
cavity temperature to 100 mK with a dilution 
refrigerator [3J. 

5. Quantum nondemolition measurements of the 
cyclotron and spin energy level are realized using 
a one-particle self-excited oscillator [llj . 
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6. The stored electron serves as its own magne- 
tometer, allowing the accumulation of lineshape 
statistics over days, revealing that a broadening 
of the expected lineshapes is the major remaining 
uncertainty [2]. 

Following sections will discuss the quantum cyclotron 

and the new methods. 



C. A Long History 

As befits one of the few properties of the electron that 
can be accurately measured, the Harvard magnetic mo- 
ment measurements detailed here are only the latest in 
a long history of measurements that make use of dif- 
ferent methods. The early history [T2] established that 
g/2 ~ 1. A series of more precise measurements followed 
at the University of Michigan, by measuring the differ- 
ence of the cyclotron and spin precession frequencies of 
keV electrons traveling on helical orbits in a magnetic 
field, concluding with a 3500 ppt measurement of g/2 
[12*] . (Here ppt refers to 1 part in 10 12 , and ppb refers to 1 
part in 10 9 .) Research groups at the University of Mainz 
and the University of Washington (UW) next developed 
methods to measure the electron magnetic moment using 
a large number of electrons stored in a Penning trap (131 — 
fl6] . Out of these efforts came the capability to suspend 
and detect a single electron in a Penning trap [T7] at the 
UW. A few years later a measurement was made with one 
electron [18]. Over the next decade these methods were 
refined, culminating in the celebrated 1987 measurement 
already mentioned [3] that reported an uncertainty of 4 
ppt. 



D. The Fine Structure Constant 



The fine structure constant, 



47reo^c' 



(3) 



gives the strength of the electromagnetic coupling in the 
low-energy limit. The energy scales for atoms are set by 
powers of a times the electron rest energy, mc 2 . For hy- 
drogen atoms the binding energy scale is a 2 mc 2 , the fine 
structure splitting scale is a 4 mc 2 , and Lamb shift scale is 
a 5 mc 2 . The quantum Hall conductance is proportional 
to the fine structure constant. The fine structure con- 
stant is already a crucial ingredient in our system of fun- 
damental constants [19j [20] and it will acquire a more 
prominent role if plans to redefine the SI system of units 
[2"T] go forward. 

Sec. IVIIBI shows how the new measurements of the 
electron g/2 determine a to be [2] [22] 123] . 



aT 1 = 137.035 999 084(51) 



[0.37 ppb]. (4) 



The uncertainty in a is now limited a bit more by the 
need for a higher-order QED calculation (underway [23]) 



than by the measurement uncertainty in g/2. The total 
0.37 ppb uncertainty in a is now about 12 times smaller 
than that of the next-most-precise independent method 
(Fig.§. 
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II. ONE-ELECTRON QUANTUM CYCLOTRON 
A. Electron in a Magnetic Field 

For an electron in a magnetic field, g/2 is specified by 
its spin and cyclotron frequencies, v 8 and v c , 
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or equivalently by their difference (the anomaly frequency 
v a = v s — v c ) and v c . Because u s and v c differ by only a 
part-per-thousand, measuring v a and v c to a precision of 
1 part in 10 10 gives g/2 to 1 part in 10 13 . 

Although we cannot measure accurately with one elec- 
tron in free space because the electron would not stay in 
one place long enough, two features of determining g/2 
are already apparent in Eq. [5] First, one can determine 
g/2 by measuring a ratio of frequencies. This is fortunate 
because there is nothing in physics that can be measured 
more accurately than a frequency (the art of time keep- 
ing being developed being so highly developed) except 
for a ratio of frequencies. Second, although both of these 
frequencies depend upon the magnetic field, the field de- 
pendence drops out of the ratio. The magnetic field thus 
needs to be stable only on the time scale on which both 
frequencies can be measured, and no absolute calibration 
of the magnetic field is required. 



B. Electron in a Penning Trap 

An ideal Penning trap confines an electron using 
a magnetic field Bz with an additional electrostatic 
quadrupole potential V ~ z 2 — p 2 /2 [28]. This poten- 
tial confines the electron axially with frequency v z and 
shifts the cyclotron frequency from the free-space value 
v c to v c . The latter frequency is also slightly shifted 
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FIG. 3. Lowest cyclotron and spin levels of an electron in a 
Penning trap. 



by the unavoidable leading imperfections of a real lab- 
oratory trap - a misalignment of the symmetry axis of 
the electrostatic quadrupole and the magnetic field, and 
quadratic distortions of the electrostatic potential. 

The lowest cyclotron energy levels (with quantum 
numbers n — 0, 1, . . .) and the spin energy levels (with 
quantum numbers m s — ±1/2) (Fig. |3| are given by 

E(n,m s ) — —hv c rn s -\-(n+^)hv c — ^hS(n+^+m s ) 2 . (6) 

The third term in Eq. [6] is the leading relativistic correc- 
tion [33] to the energy levels. Special relativity makes 
the transition frequency between two cyclotron levels 
\n, m s ) <-> \n + 1, m s ) decrease from v c to D c + Az/ C , with 
the shift 



circular orbit for which the electric field for the trap and 
the motional magnetic field exactly cancel. 

The invariance theorem applies for a perfect Penning 
trap, but also in the presence of the mentioned imperfec- 
tion shifts of the eigenfrequencies for a real trap. This 
theorem, together with the well-defined hierarchy of trap 
eigenfrequencies, v c ^> v z 3> v m ^> 8, yields an approx- 
imate expression that is sufficient at our accuracy. We 
thus determine the electron g/2 using 



Da - 9 2 J(2f c ) 
U + 3<5/2 + ?2/(2/c) 



(10) 



The determination requires four inputs. First and sec- 
ond are high-precision measurements of the transition 
frequencies 
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represented by the red and blue arrows in Fig. [6j Third 
is a relatively lower precision measurement of the axial 
frequency v z . Fourth is the cavity shift Ag cav /2 that 
arises from the interaction of the cyclotron motion and 
the trap cavity and is discussed in detail in Sec. \V\ 
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depending upon the spin state and cyclotron state. This 
very small shift, with 



5/v c — hv c / (mc 2 ) w 10 



(8) 



is nonetheless significant at our precision. However, an 
essentially exact treatment of the relativsitic shift is pos- 
sible because single quantum transitions are resolved. 
The relativistic shift thus contributes no uncertainty to 
our measurement. This is a key advantage of the quan- 
tum cyclotron over previous measurements systems [1], 
in which an unknown distribution of cyclotron states was 
excited each with a different relativistic shift. 

To determine g/2, we must rewrite Eq. [5] in terms 
of measurable frequencies of an electron bound in the 
trap. The needed free-space cyclotron frequency, v c — 
eB/{2%m), is deduced by use of the Brown-Gabrielse in- 
variance theorem 1301. 



(^c) 2 = (^) 2 + (^) 2 + (P m ) i 
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The three measurable eigenfrequencies on the right in- 
clude the cyclotron frequency v c for the quantum cy- 
clotron motion we have been discussing. The second 
measurable eigenfrequency is the axial oscillation fre- 
quency v z for the nearly-harmonic, classical electron mo- 
tion along the direction of the magnetic field. The third 
measurable eigenfrequency is the magnetron oscillation 
frequency for the classical magnetron motion along the 



III. EXPERIMENTAL REALIZATION 

A. Cylindrical Penning Trap 

A cylindrical Penning trap (Fig. |4| is the key device 
that makes these measurements possible. It was invented 
[6] and demonstrated [7] to provide boundary conditions 
that produce a controllable and understandable radia- 
tion field within the trap cavity. Spontaneous emission 
can be significantly inhibited at the same time as cor- 
responding shifts of the electron's oscillation frequencies 
are avoided. The latter has not been possible [3T] with 
the hyperbolic Penning traps of earlier experiments [1], 
which have electrodes approximating the equipotentials 
of an electrostatic quadrupole. 
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FIG. 4. Cylindrical Penning trap cavity used to confine a 
single electron and inhibit spontaneous emission. 

The first function of the trap electrodes is to produce a 
very good approximation to an electrostatic quadrupole 
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potential. This requires careful choice of the relative ge- 
ometry of the electrodes [6] . The electrodes of the cylin- 
drical trap are symmetric under rotations about the cen- 
ter axis (z), which is parallel to the spatially uniform 
magnetic field (Bz). The potential (about 100 V) ap- 
plied between the endcap electrodes and the ring elec- 
trode provides the basic trapping potential and sets the 
axial frequency v z of the nearly harmonic oscillation of 
the electron parallel to the magnetic field. The poten- 
tial applied to the compensation electrodes is adjusted 
to tune the shape of the potential, to make the oscilla- 
tion as harmonic as possible. The tuning does not change 
v z very much owing to an orthogonalization [6., J2] that 
arises from the geometry choice. What we found was 
that one electron could be observed within a cylindrical 
Penning trap with as good or better signal-to-noise ratio 
than was realized in hyperbolic Penning traps. 



that are indexed with integers 



m = 0,l,2,--- 
n=l,2,.-- 
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and are functions of the nth zeros of Bessel functions and 
their derivatives 
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The zeros force the boundary conditions at the cylindrical 
wall. All but the m = modes are doubly degenerate. 

Of primary concern is the magnitude of the transverse 
electric fields since only these components couple to cy- 
clotron motion. For both TE and TM modes, the trans- 
verse components of E are proportional to 



B. Radiation Field in a Cylindrical Trap 

The second function of the trap electrodes is to form 
a microwave cavity whose radiation properties are well 
understood and controlled. The density of states for a 
perfect right circular cylinder are the familiar transverse 
electric (TE) and transverse magnetic (TM) radiation 
modes for such a geometry. In the real trap cavity, the 
perturbation caused by the small space between the elec- 
trodes is minimized by the use of "choke flanges" - small 
channels that tend to reflect the radiation leaking out 
of the trap back to cancel itself, and thus to minimize 
the losses from the trap. The measured radiation modes 
(Fig. 10) are close enough to the calculated frequencies 



for a perfect cylindrical cavity that we have been able 
to identify more than 100 different radiation modes for 
such trap cavities [8]. The spatial properties of the elec- 
tric and magnetic field for the radiation that builds up 
within the cavity are thus quite well understood. Some 
of the modes couple to cyclotron motion of an electron 
centered in the cavity, others couple to the spin of a cen- 
tered electron, and still others have the symmetry that 
we hope will one day allow us to sideband cool the axial 
motion. 

A right circular cylinder of diameter 2/?o and height 
2zq admits two classes of electromagnetic fields, trans- 
verse electric TE„ mp and transverse magnetic TM mnp , 
each identified by three indices that describe their spe- 
cific geometry (see e.g. [Ml Sec. 8.7]). These TE and TM 
modes have characteristic frequencies, 
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for even p, 

(.^(p-iya cos(gl) for odd p. 

(17) 

For an electron close to the cavity center, (z 0), only 
modes with odd p thus have any appreciable coupling. 

The transverse components of the electric fields are 
also proportional to either the order-m Bessel function 
times m/p or to the derivative of the order-m Bessel func- 
tion. Close to the cavity center (p ss 0), 
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In the limit p — > 0, all but the m = 1 modes vanish. 

For a perfect cylindrical cavity the only radiation 
modes that couple to an electron perfectly centered in 
the cavity are TE ln ( odd ) and TM ln ( odd ). If the electron 
is moved slightly off center axially it will begin to couple 
to radiation modes with mnp = ln(even). If the electron 
is moved slightly off-center radially it similarly begins to 
couple to modes with ra/1. 
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100 mK and 5 T 



The trap cavity is cooled to 0.1 K or below via thermal 
contact with the mixing chamber of an Oxford Instru- 
ments Kelvinox 300 dilution refrigerator (Fig. [5| . They 
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FIG. 5. The apparatus. The solenoid and electrodes that 
form the Penning trap are in red. The dilution refrigerator is 
in blue. Cryogen spaces are hatched. 



are housed within a separate vacuum enclosure that is en- 
tirely at the base temperature. Measurements on an ap- 
paratus with a similar design but at 4.2 K found the vac- 
uum in the enclosure to be better than 5 x 10 -17 torr [54] . 
Our much lower temperature should make our back- 
ground gas pressure much lower. We are able to keep 
one electron suspended in our apparatus for as long as 
desired - regularly months at a time. Substantial reser- 
voirs for liquid helium and liquid nitrogen make it possi- 
ble to keep the trap cold for five to seven days before the 
disruption of adding more cryogens is required. 

The trap electrodes and their vacuum container are 
located within a superconducting solenoid (Fig. [5]) that 
makes a very homogeneous magnetic field over the in- 
terior volume of the trap cavity. A large dewar sitting 
on top of the solenoid dewar provides the helium needed 
around the dilution refrigerator below. The supercon- 
ducting solenoid is entirely self-contained, with a bore 
that can operate from room temperature down to 77 K. 
It possesses shim coils capable of creating a field homo- 
geneity better than a part in 10 s over a 1 cm diameter 
sphere and has a passive "shield" coil that reduces fluc- 
tuations in the ambient magnetic field [35] ES]- When 
properly energized (and after the steps described in the 
next section have been taken) it achieves field stability 
better than a part in 10 9 per hour. We regularly observe 
drifts below 10 -9 per night. 



D. Attaining High Stability 

Measuring the electron g/2 with a precision of parts in 
10 13 requires careful attention to making a stable trap- 



ping potential. Even more important is a stable magnetic 
field since the frequencies f c and v a that we measure are 
both proportional to B, and we are not able to measure 
these frequencies at exactly the same time. 

A major defense against external field fluctuations is a 
high magnetic field. This makes fluctuations from outside 
sources relatively smaller. The largest source of ambient 
magnetic noise is a subway that produces 50 nT (0.5 
mG, 10 ppb) fluctuations in our lab and that would limit 
us to four hours of data taking per day (when the sub- 
way stops running) if we did not shield the electron from 
them. Eddy currents in the high-conductivity aluminum 
and copper cylinders of the dewars and the magnet bore 
shield high-frequency fluctuations [37]. For slower fluc- 
tuations, the aforementioned shelf-shielding solenoid [35] 
has the correct geometry to make the central field always 
equal to the average field over the solenoid cross-section. 
This translates flux conservation into central-field con- 
servation, shielding external fluctuations by more than a 
factor of 150 [35]. 

Stabilizing the field produced by the solenoid requires 
that care is taken when the field value is changed, since 
changing the current in the solenoid alters the forces 
between windings. Resulting stresses can take months 
to stabilize if the coil is not pre-stressed by "over- 
currenting" the magnet. Our recipe is to overshoot the 
target value by a few percent of the change, undershoot 
by a similar amount, and then move to the desired field, 
pausing several minutes after each change. 

The apparatus in Fig.[5jevolved historically rather than 
being designed for maximum magnetic field stability in 
the final configuration. Because the solenoid and the trap 
electrodes are suspended from widely separated support 
points, temperature and pressure changes can cause the 
electrodes to move relative to the solenoid. Apparatus 
vibrations can do the same. Insofar as the magnetic field 
is not perfectly homogeneous, despite careful adjusting 
of the persistent currents in ten superconducting shim 
coils, such relative motion changes the field seen by the 
electron. 

To counteract this, we have long regulated the five He 



TABLE I. Typical trap parameters as well as frequencies and 
damping rates for each degree of freedom. The damping rates 
include coupling to a detection circuit for "/ z and inhibited 
spontaneous emission for "/ c . To sample the radiation modes 
of the electrode cavity, we change B, and hence u c , -y c , and 
v s (see Sec. [Vf. 
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and N2 pressures in the cryostats to 50 ppm to main- 
tain the temperature of both the bath and the solenoid 
itself [38] EH] • Recently we also relocated the dilution re- 
frigerator vacuum pumps to an isolated room at the end 
of a 12 m pipe run. This reduced vibration by more than 
an order of magnitude at frequencies related to the pump 
motion and reduced the noise level for the experimenters 
but did not obviously improve the g/2 data. 

Because some of the structure establishing the rela- 
tive location of the trap electrodes and the solenoid is 
at room temperature, changes in room temperature can 
move the electron in the magnetic field. The lab tem- 
perature routinely cycles 1 - 2 K daily, so we house the 
apparatus in a large, insulated enclosure within which we 
actively regulate the air temperature to 0.1 K. A refrig- 
erated circulating bath (ThermoNeslab RTE-17) pumps 
water into the regulated zone and through an automo- 
bile transmission fluid radiator, heating and cooling the 
water to maintain constant air temperature. Fans couple 
the water and air temperatures and keep a uniform air 
temperature throughout. 

The choice of materials for the trap electrodes and its 
vacuum container is also crucial to attaining high field 
stability [TJ 2D]. Copper trap electrodes, for example, 
have a nuclear paramagnetism at 0.1 mK that makes the 
electron see a magnetic field that changes at an unac- 
ceptable level with very small changes in trap tempera- 
ture. We thus use only low-Curie-constant materials such 
as silver, quartz, titanium, and molybdenum at the re- 
frigerator base temperature and we regulate the mixing 
chamber temperature to better than 1 mK. 

A stable axial frequency is also extremely important 
since small changes in the measured axial frequency re- 
veal one-quantum transitions of the cyclotron and spin 



energy (as will be discussed in Sec. IIIF). A trapping 



potential without thermal fluctuations is provided by a 
charged capacitor (10 /iF) that has a very low leakage re- 
sistance at low temperature. We add to or subtract from 
the charge on the capacitor using 50 ms current pulses 
sent to the capacitor through a 100 MSI resistor as needed 
to keep the measured axial frequency constant. Because 
of the orthogonalized trap design [B] already discussed, 
the potential applied to the compensation electrodes (to 
make the electron see as close to a pure electrostatic 
quadrupole potential as possible) has little effect upon 
the axial frequency. 



E. One Electron: Its Motions and Damping 

We load a single electron using an electron beam from 
an atomically sharp tungsten field-emission tip. A hole 
in the bottom endcap electrode admits the beam, which 
hits the top endcap electrode and releases gas atoms cry- 
opumped on the surface. Collisions between the beam 
and gas atoms allows an electron to fall into the trap. 
Adjusting the beam energy and the time it is left on con- 
trols the number of electrons loaded. 



The electron has three motions in the Penning trap 
formed by the B — 5.4 T magnetic field, and the electro- 
static quadrupole potential. The cyclotron motion in the 
trap has a frequency D c ss 150 GHz. The axial frequency, 
for the harmonic oscillator parallel to the magnetic field 
direction, is v z « 200 MHz. A circular magnetron mo- 
tion, perpendicular to B, has an oscillation frequency, 
v m ~ 133 kHz. The spin precession frequency, which we 
do not measure directly, is a part-per-thousand higher 
than the cyclotron frequency. The frequency difference 
is the anomaly frequency, v a ~ 174 MHz, which we do 
measure directly. 

The undamped spin motion is essentially uncoupled 
from its environment |28j . The cyclotron motion is only 
weakly damped. By controlling the cyclotron frequency 
relative to that of the cavity radiation modes, we alter 
the density of radiation states and inhibit the sponta- 
neous emission of synchrotron radiation [551 ODE] by 10 to 
50 times the (90 ms) -1 free-space rate. Blackbody pho- 
tons that could excite from the cyclotron ground state 
arc eliminated because the trap cavity is cooled by the 
dilution refrigerator to 100 mK [3J. The axial motion is 
cooled by a resonant circuit at a rate j z ~ (0.2 s) _1 to 
as low as 230 mK (from 5 K) when the detection ampli- 
fier is off. The magnetron radius is minimized with axial 
sideband cooling |28j . 



F. QND Detection 

Quantum nondemolition (QND) detection has the 
property that repeated measurements of the energy 
eigenstate of the quantum system will not change the 
state [42] . This is crucial for our detection of one- 
quantum transitions in the cyclotron motion where we 
do not want the detection system to be producing the 
transitions that we observe. In this section we discuss 
the QND coupling and in the next our readout system. 

Detecting a single 150 GHz photon from the decay of 
one cyclotron energy level to the level below is very dif- 
ficult - because the frequency is so high and because it 
is difficult to cover the solid angle into which the pho- 
ton could be emitted. Instead we get the single-photon 
sensitivity by coupling the cyclotron motion to the or- 
thogonal axial motion at 200 MHz, a frequency at which 
we are able to make sensitive detection electronics |43j . 
The QND nature of the detection means the thermally 
driven axial motion of the electron does not change the 
state of the cyclotron motion. 

We use a magnetic bottle gradient that is familiar from 
plasma physics and from earlier electron measurements 

I!, 



AB = B 2 [(i 



72) 



zpp 



(19) 



with B2 = 1540 T/m 2 . The gradient arises from a pair 
of thin nickel rings (Fig. [4]) that are completely saturated 
in the strong field from the superconducting solenoid. To 
lowest order the rings modify B by ss —0.7% - merely 
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changing the magnetic field that the electron experiences 
without affecting our measurement. 

The formal requirement for a QND measurement is 
that the Hamiltonian of the quantum system (i.e. the cy- 
clotron Hamiltonian) and the Hamiltonian describing the 
interaction of the quantum system and the classical mea- 
surement system must commute. The Hamiltonian that 
couples the quantum cyclotron and spin motions to the 
axial motion does so. It has the form — /j,B, where \x is the 
magnetic moment associated with the cyclotron motion 
or the spin. The coupling Hamiltonian thus has a term 
that goes as fiz 2 . This term has the same spatial sym- 
metry as does the axial Hamiltonian, H = ^m(2ni' z ) 2 z 2 . 
A change in the magnetic moment that takes place from 
a one-quantum change in the cyclotron or spin magnetic 
moment thus changes the observed axial frequency of the 
suspended electron. 
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FIG. 6. Two quantum leaps: A cyclotron jump (a) and spin 
flip (b) measured in a QND manner through shifts in the axial 
frequency. 

The result is that the frequency of the axial motion v z 
shifts by 



= S B (n + m s ), 



(20) 



in proportion to the cyclotron quantum number n and 
the spin quantum number m s . Fig. [6] shows the AD Z — 4 
Hz shift in the 200 MHz axial frequency that takes place 
for one-quantum changes in cyclotron (Fig. [6^,) and spin 
energy (Fig. [6Jd). The 20 ppb shift is easy to observed 
with an averaging time of only 0.5 s. We typically mea- 
sure with an averaging time that is half this value. 



G. One-Electron Self-Excited Oscillator 

Cyclotron excitations and spin flips are generally in- 
duced while the detection system is off, as will be dis- 
cussed. After an attempt to excite the cyclotron motion 
or to flip the spin has been made, the detection system is 
then turned on to detect the state of the system. Spon- 
taneous emission of synchrotron radiation from the cy- 



clotron motion is inhibited (Sec. IIIH) to give the time 
that is needed. The small axial frequency shifts which 
signal changes in cyclotron quantum number or spin are 
measured at the required precision using a one-electron 
self-excited oscillator [11] that turns on rapidly. 

The 200 MHz axial frequency lies in the radio- 
frequency (rf) range which is more experimentally ac- 
cessible than the microwave range of the 150 GHz cy- 



clotron and spin frequencies, as mentioned. Neverthe- 
less, standard rf techniques must be carefully tailored for 
our low-noise, cryogenic experiment. The electron axial 
oscillation induces image currents in the trap electrodes 
that are proportional to the axial velocity of the elec- 
tron [28l |45] . An inductor is placed in parallel with the 
capacitance between two trap electrodes to cancel the re- 
actance of the capacitor which would otherwise short out 
the induced signal. The rf loss in the tuned circuit that 
is formed is an effective resistance that damps the axial 
motion. 

The voltage that the electron motion induces across 
this effective resistance is amplified with two cryogenic 
amplifiers. The heart of each amplifier is a single-gate 
high electron mobility transistor (Fujitsu FHX13LG). 

The first amplifier is at the 100 mK dilution refriger- 
ator base temperature. Operating this amplifier with- 
out crashing the dilution refrigerator requires operating 
with a power dissipation in the FET that is three or- 
ders of magnitude below the transistor's 10 mW design 
dissipation. The effective axial temperature for the elec- 
tron while current is flowing through the FET is about 
5 K, well above the ambient temperature. Very care- 
ful heat sinking makes it possible for the effective axial 
temperature of the electron to cool to below 350 mK in 
approximately one second, taking the electron axial mo- 
tion to this temperature. Cyclotron excitations and spin 
flips are induced only when the axial motion is so cooled, 
since the electron is then making the smallest possible 
excursion in the magnetic field gradient. 

The second amplifier is mounted on the nominally 600 
mK still of the dilution refrigerator. This amplifier coun- 
teracts the attenuation of thermally-isolating but lossy 
stainless steel transmission line that caries the ampli- 
fied signal out of the refrigerator. The second amplifier 
boosts the signal above the noise floor of the first room- 
temperature amplifier. 

We feed this signal back to trap electrodes as a drive. 
Because the induced signal is proportional to the elec- 
tron's axial velocity, this feedback alters the axial damp- 
ing force, a force that is also proportional to the elec- 
tron velocity. Changing the feedback gain thus changes 
the damping rate. As the gain increases, the damping 
rate decreases as does the effective axial temperature of 
the electron, in accord with the fluctuation dissipation 
theorem [46] . The invariant ratio of the separately mea- 
sured damping rate and the effective temperature has 
been demonstrated [UJ , thereby also demonstrating that 
the amplifier adds very little noise to the feedback. 

Setting the feedback gain to make the drive exactly 
cancel the damping in the attached circuit could sustain 
a large axial oscillation amplitude, in principle. However, 
the gain cannot be perfectly adjusted and noise fluctua- 
tions will always drive the axial oscillation exponentially 
away from equilibrium. We thus stabilize the oscillation 
amplitude using a digital signal processor (DSP) that 
Fourier transforms the signal, and adjusts the feedback 
gain in real time to keep the signal at a fixed value. 
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This one-particle self-excited oscillator is turned on af- 
ter an attempt has been made to excite the cyclotron 
energy up one level, or to flip the spin. The frequency 
of the axial oscillation that rapidly stabilizes at a large 
and easily detected amplitude is then measured. Small 
shifts in this frequency reveal whether the cyclotron mo- 
tion has been excited or whether the spin has flipped, as 
illustrated in Fig. [6] 



H. Inhibited Spontaneous Emission 

One of the early papers in what has come to be know 
as cavity electrodynamics was an observation of inhib- 
ited spontaneous emission within a Penning trap |41j - 
the first time that inhibited spontaneous emission was 
observed within a cavity and with only one particle - 
as anticipated earlier [J5J As already mentioned, 

the cylindrical Penning trap [6] was invented to provide 
boundary conditions that would allow the control of the 
electron-cavity coupling, using an understandable geom- 
etry that allows the calculation of cavity shifts to the 
electron's cyclotron frequency. 

The spontaneous emission rate can be easily measured 
directly, by making a histogram of the time the electron 
spends in the first excited state after being excited by 
a microwave drive injected into the trap cavity. Fig. [JJ 
shows a sample histogram which fits well to an exponen- 
tial (solid curve) with a lifetime of 3.4 s in this example. 



The spontaneous emission rate into free space is 



1 4e 2 
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FIG. 7. Histogram of the time that the electron spends in the 
first excited state with an exponential fit. The decay time, 
3.4(1) s in this example, depends on how close the cyclotron 
frequency is to neighboring radiation modes of the trap cavity. 
Lifetimes as long as 16 s have been observed. 

Stimulated emission is avoided by making these obser- 
vations only when the cavity is at low temperature so 
that effectively no black body photons are present. The 
detector makes thermal fluctuations of the axial oscil- 
lation amplitude, and these in turn make the cyclotron 
frequency fluctuate. For measuring the cyclotron decay 
time, however, this does not matter because the fluctua- 
tions in axial amplitude are small compared to the 2 mm 
wavelength of the radiation that excites the cyclotron 
motion. 



47reo 3mc 3 i 



ms 



(21) 



The measured rate in the example of Fig. [7] is thus sup- 
pressed by a factor of 38. The density of states within the 
cylindrical trap cavity is not that of free space. Instead 
the density of states for the radiation is peaked at the res- 
onance frequencies of the radiation modes of the cavity, 
and falls to very low values between the radiation modes. 
We attain the inhibited spontaneous emission by tuning 
the magnetic field so that the cyclotron frequency is far 
from radiation modes. With the right choice of magnetic 
field we have increased the lifetime to 16 s, which is a 
cavity suppression of spontaneous emission by a factor of 
180. 



In Sec. VD we report on using the direct measure- 
ments of the radiation rate for electron cyclotron motion 
to probe the radiation modes of the cavity, with the radi- 
ation rate increasing sharply at frequencies that approach 
a resonant mode of the cavity. 



I. Quantum Jump Spectroscopy 

We probe the cyclotron and anomaly resonances with 
quantum jump spectroscopy, in which we apply a drive in 
discrete frequency steps, checking between applications 
for a one-quantum transition and building a histogram 
of the ratio of excitations to attempts at each frequency. 
Fig. [9] show the observed lineshapes upon which our best 
measurement is based. We discuss how the measured 
points are obtained first, and then discuss theoretical 
lineshapes in Sec. |IV| 

A typical data run consists of alternating scans of the 
cyclotron and anomaly lines and occurs at night, with 
daytime runs only possible on Sundays and holidays when 
the ambient magnetic field noise is lower. Interleaved ev- 
ery three hours among these scans are periods of mag- 
netic field monitoring to track long-term drifts using the 
electron itself as the magnetometer. In addition, we con- 
tinuously monitor over fifty environmental parameters 
such as refrigerator temperatures, cryogen pressures and 
flows, and the ambient magnetic field in the lab so that 
we may screen data for abnormal conditions and trou- 
bleshoot problems. 

Cyclotron transitions are driven by injecting mi- 
crowaves into the cavity. The microwaves originate as a 
15 GHz drive from a signal generator (Agilent E8251A) 
whose low-phase-noise, 10 MHz oven-controlled crys- 
tal oscillator serves as the timebase for all frequencies 
in the experiment. After passing through a waveg- 
uide that removes all subharmonics, the signal enters 
a microwave circuit that includes an impact ionization 
avalanche transit-time (IMPATT) diode, which multi- 
plies the frequency by ten and outputs the f c drive at 
a power of 2 mW. Voltage-controlled attenuators reduce 
the strength of the drive, which is broadcast from a 
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room temperature horn through teflon lenses to a horn 
at 100 mK (Fig. [5]) and enters the trap cavity through an 
inlet waveguide (Fig. H). 

Anomaly transitions are driven by potentials, oscillat- 
ing near v ai applied to electrodes to drive off-resonant ax- 
ial motion through the magnetic bottle gradient (Eq. 19 ). 
The gradient's zpp term mixes the driven oscillation of 
z at v a with that of p at f c to produce an oscillating 
magnetic field perpendicular to B as needed to flip the 
spin. The axial amplitude required to produce the de- 
sired transition probability is too small to affect the line- 



shape (Sec. VI); nevertheless, we apply a detuned drive 



of the same strength during cyclotron attempts so the 
electron samples the same magnetic gradient. 

Quantum jump spectroscopy of each resonance follows 
the same procedure. With the electron prepared in the 
spin-up ground state |0, |), the magnetron radius is re- 
duced with 1.5 s of strong sideband cooling at v z + v m 
with the SEO turned off immediately and the detection 
amplifiers turned off after 0.5 s. After an additional 1 s 
to allow the axial motion to thermalize with the tuned 
circuit, we apply a 2 s pulse of either a cyclotron drive 
near f c or an anomaly drive near T> a with the other drive 
applied simultaneously but detuned far from resonance. 
The detection electronics and SEO are turned back on; 
after waiting 1 s to build a steady-state axial amplitude, 
we measure v z and look for a 20 ppb shift up (from a cy- 
clotron transition) or down (from an anomaly transition 
followed by a spontaneous decay to |0, — |)) in frequency. 
Cavity-inhibited spontaneous emission provides the time 
needed to observe cyclotron transitions before decay. The 
several-cyclotron-lifetimes wait for a spontaneous decay 
after an anomaly attempt is the rate-limiting step in the 
spectroscopy. After a successful anomaly transition and 
decay, simultaneous cyclotron and anomaly drives pump 
the electron back to |0, |). All timing is done in hard- 
ware with a pulse generator. We probe each resonance 
line with discrete excitation attempts spaced in frequency 
by approximately 10% of the linewidth. We step through 
each drive frequency on the f c line, then each on the D a 
line, and repeat. 



lows the normalization of the raw cyclotron and anomaly 
data to a common magnetic field. 



K. Measuring the axial frequency 



Our expression for g/2 (Eq. |10|) requires a measure- 
ment of the axial frequency v z . For a relative uncer- 
tainty in g below 0.1 ppt, we must know v z to better 
than 50 ppb, or 10 Hz. This is easily done. 

While we routinely measure v z during QND detection 
of cyclotron and spin states, this self-excited oscillation 
frequency includes an amplitude-dependent anharmonic 
shift from the low-amplitude, thermally excited axial mo- 
tion during cyclotron and anomaly excitation attempts. 
The shift is typically a few hertz. We account for this 
shift by directly measuring the thermal axial frequency 
with the amplifiers on; the axial resonance appears as a 
narrow dip where the electron has "shorted-out" the am- 
plifier noise |45j . This dip frequency is negligibly differ- 
ent than the slightly lower-amplitude one that pertains 
during cyclotron and anomaly excitation, which occurs 
with the amplifiers off. All additional shifts from inter- 
action with the amplifier or anomaly-drive- induced power 
shifts [150] are negligible at our precision. Our result for 
g/2 has no uncertainty arising from the measurement of 



IV. CYCLOTRON AND ANOMALY 
LINESHAPES 

Quantum jump spectroscopy resolves the cyclotron 
and anomaly resonance lines, and we rely on an invari- 
ant property of the expected lineshapes (their mean fre- 
quency) to determine f c and v a with uncertainties less 
than the linewidths. The observed lines are slightly 
broader than expected, an effect we attribute to magnetic 
field fluctuations. Our understanding of this broadening 
is the primary limitation of the most accurate of our mea- 
surements. 



J. The Electron as Magnetometer 

Slow drifts of the magnetic field are corrected using 
the electron itself as a magnetometer. Accounting for 
these drifts allows the combination of data taken over 
many days, giving a lineshape signal-to-noise that allows 
the systematic investigation of lineshape uncertainties at 
each field. For a half-hour at the beginning and end of 
a run and again every three hours throughout, we alter 
our cyclotron spectroscopy routine by applying a stronger 
drive at a frequency below f c . Using the same timing as 
above but a ten-times-finer frequency step, we increase 
the drive frequency until observing a successful transi- 
tion. We then jump back 60 steps and begin again. Fit- 
ting a polynomial to such cyclotron- edge tracking data al- 



A. Lineshape model 

Detailed discussion of the cyclotron and anomaly res- 
onance lineshapes are available [28, 43j[5T]. Here we give 
a summary, stressing what must be assumed and calcu- 
lated to enable the measurement of the electron magnetic 
moment and the uncertainties in the measurement. 

For cyclotron and anomaly drives left on for a time 
much longer than the inverse-linewidth and inverse-axial- 
damping rate 7~ , the probability P for a transition to 
occur after a time T is [5T] 

P= i{l-exp[-7r™ 2 X H]}, (22) 
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which depends on the Rabi frequency Q and the lineshape 
function x( w ) an< i saturates at 1/2 for strong drives. The 
radiative decay of cyclotron excitations reduces the cy- 
clotron lineshape saturation value. For a delay t between 
the end of the drive and the beginning of the state mea- 
surement, the lineshape saturates at exp(— j c t)/2. 

In general, the lineshape x( w ) is the Fourier transform 
of a correlation function which is related to the 

statistical average of any fluctuations in the magnetic 
field, uj(t) [ST]: 



exp 



dt'u)(t') 



(23) 
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1. During quantum jump spectroscopy 
motion 



Brownian axial 



FIG. 8. The theoretical Brownian-motion lineshape for vari- 
ous j z /Auj and 7c/7z- Our cyclotron line has 7 z /Acj ~ 10 -2 . 
Our anomaly line has 7 z /Ao; ~ 10. For our j z , 7 c /7z = 10 _1 - 
10 -2 correspond to lifetimes of 1.6 s - 16 s. The infinite- 
lifetime limit is 7c/7z = 0. 



The same magnetic bottle that makes a QND coupling 
of the cyclotron and spin energies to the axial frequency 
couples the axial energy to the cyclotron and anomaly 
frequencies and is the primary source of the observed 
lineshape. For an electron on axis (p = 0), it adds a 
z 2 dependence to the magnetic field and thus to the cy- 
clotron and anomaly frequencies, here collectively w: 



uj(z) = UJq 1 + 



B 
~B 



(24) 



When probing the cyclotron or anomaly resonance for 
quantum jump spectroscopy, the axial motion is in ther- 
mal equilibrium with the detection amplifier. Thus, the 
axial position z undergoes Brownian motion. The line- 
shape x( w ) is a statistical average of these Brownian fluc- 
tuations and takes various forms depending on the rel- 
ative lengths of the fluctuation timescale — the inverse- 
axial-damping time 7J 1 — and the inverse- linewidth co- 
herence time. This linewidth roughly corresponds to the 
frequency shift at the root-mean-square thermal axial 
amplitude: 



B 2 2 B 2 kT z 

Aw = uj — z lms = uj — - . 

B B mujt 



(25) 



The relevant correlation function x(t) is found by in- 
serting Eq. [24] whose time-dependence is in the axial 
Brownian motion, into Eq. [23] to get 



lutions for the lineshape [51] . 



XM = -Re 

7T 



where 



7 lz 



i(w-w )tg-i(7' — u)tp-^1ct 



e 2 



< /o ^(7' + 7,) 2 -(7'-7,) 2 e-^ 



(27a) 



-Re 



7T 



7 7z 



(7' + 7,) 2 

(7'-7,) 2 "(7'+7,r 2rl 



= --Re 

7T 



^( 7z + 7 ') 2 



x 2 F 1 [1,-K;1- K: 



ilz - if 

{lz+i? 



K = 



7' = + 4i7 2 Aw, 

2i(iO - UJq) + ~1z-j -1c 

27' 



(27b) 

(27c) 

(28) 
(29) 



"Re" denotes the real part, and 2-Fi(a, b; c; z) is a hyper- 
geometric function. Examples of the lineshape for various 
values of 7 2 /Aw and j z /lc are shown in Fig. [8] 

In contrast to some previous presentations of the line- 
shape [55] HSJi we have not taken the limit of low cy- 
clotron damping^] Our resolution of the anomaly line is 



X(t) = e'^ 



exp 



(26) 



Taking the statistical average gives three equivalent so- 



1 Brown does keep a non-zero cyclotron damping rate when deriv- 
ing the anomaly line, with a result 1511 Eq. 6.12] identical to that 
presented here. The derivation assumes an anomaly excitation 
technique that differs from ours (Sec. |III but the only effect is 
a redefinition of the Rabi frequency 50 . 
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fine enough that we must include the broadening from 
the finite lifetime of the |l, — 5) state. 

The cyclotron and anomaly lines are in two limits of 
this lineshape. For the cyclotron line, 7 Z /Aw « 10~ 2 <C 
1 and the axial motion is essentially decoupled from the 
amplifier during the inverse-linewidth coherence time. 
During that time, the electron remains in a single axial 
state and the lineshape is a Lorentzian with the natu- 
ral linewidth, 7 C , and centered on the frequency given 
by Eq. [24] with the rms axial amplitude of that state as 
z. We do not know which axial state that is, however, 
and since excitation attempts occur on timescales longer 
than 7J 1 , subsequent attempts will be in different states. 
Thus, the composite lineshape after many attempts is 
the convolution of the instantaneous lineshape (the nar- 
row Lorentzian) and the Boltzmann distribution of axial 
states. That is, the lineshape is a decaying exponential 
with a sharp edge at loq and a width of Aw. The cyclotron 
line is close to this "exponential" limit and should have 
a sharp edge at the zero-axial-amplitude cyclotron fre- 
quency that is useful for quick field measurements such 
as tracking drifts. 

In contrast, the anomaly line has 7 Z /Aw ra 10 > 1, 
and the axial motion is strongly coupled to the ampli- 
fier. During the inverse-linewidth coherence time, the 
axial amplitude relaxes to the thermal z rms yielding a 
lineshape that approaches a natural-linewidth Lorentzian 
offset from loq by Aw through a Lorentzian with a width 
of 7c + 2Ao; 2 /7 2 . 



2. During cyclotron lifetime measurements - driven axial 
motion 



3. Other magnetic field fluctuations 

Although the cyclotron line is in the exponential line- 
shape limit ('jz/Aio <C 1) and should have a sharp low- 
frequency edge, all our data show an edge-width of 0.5- 
1 ppb (Fig. [9]). We model this discrepancy as addi- 
tional fluctuations in the magnetic field. (Some potential 
sources are discussed at the end of this section.) Such 
fluctuations can be described by adding a noise term rj (t) 
to the Brownian motion axial fluctuations of Eq. [24] 



w(t)=W0 [l + -j^z{tf +rj(t) 



(32) 



so that the lineshape is given by the Fourier transform of 



exp 



B 9 



1 j dfw f 1 + ~z(t'f +V(t') 



(33) 

For magnetic field noise that is not correlated with the 
axial fluctuations, the average factors into 



X{t) =e- lulot { exp 



x ( exp 



-lUJQ 



B2 

B 



dt'z(t'f 



(34) 



dt'r)(t') 



The first two factors are the Brownian-motion lineshape 
and the third is an additional noise broadening. Because 
of the Fourier transform convolution theorem, the result- 
ing noisy lineshape is the noise-free lineshape of Eq. 27 
convolved with a noise function. 



The axial motion is in thermal equilibrium during the 
resonance probes of a g/2 measurement, but it is nec- 
essarily driven during QND detection. Probing the line- 
shape with this detection drive on adds a coherent oscilla- 
tion of z to the Brownian motion considered above. The 
magnetic bottle coupling again translates this motion 
into a lineshape that depends on the amplitude of the ax- 
ial oscillation. In the weak-coupling limit (7 Z /Aw <C 1) 
that corresponds to our cyclotron line, the lineshape 
is [EQ Eq. 7.18] 



2 v /(w-w )A d w 

Iq{ — 7 - w o 

Aw 

1 w - w + A d w , 

exp(- 



(30) 



where 



Aw 



A d w = w 



_Wo_ 

Aw 



B 2 A 2 



4- Invariance of the mean frequency 

An important feature of the Brownian-motion line- 
shape, which we use in our primary line-splitting tech- 
nique, is the independence of its mean from 7 Z . For the 
low drive strengths used in this measurement, we may 
expand Eq. [22] to lowest order, 



P = ^vrra 2 x(w) 



(35) 



such that the excitation probability is linear in the line- 
shape function. Then, the average frequency of the line- 
shape, (w) = w%(w) cLo, always corresponds to that 
given by the thermal z rms in the magnetic bottle field 
(Eq.[24]). That is [EQ Eq. 1.29], 



(w) = wq + Aw. 



(36) 



B 2 



(31) 



This mean is easily verified in the exponential and 
Lorentzian limits described above, but applies for all 
7 Z /Aw. 

Importantly, any additional magnetic field noise rj 
affects both cyclotron and anomaly lines identically 
in calibrating the axial oscillation amplitude in Sec. VD (Eq. 34). If the noise fluctuates symmetrically, that is, 



A is the driven axial amplitude, 6(x) is the Heaviside 
step function, and Io(x) is the order-zero modified Bessel 
function. We make extensive use of this driven lineshape 
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the mean frequency of the noise function is zero, then the 
mean frequencies of the lines are unchanged. A non-zero 
noise function average would shift both lines by propor- 
tionally the same amount; this shift would then cancel in 
the calculation of g/2 (Eq. 10 1. 



The large-drive saturation of Eq. 22 invalidates this 
mean-frequency invariance, but we keep our excitation 
probabilities below 20%, where we expect any saturation- 
shift of the mean frequency to be smaller than the sta- 
tistical uncertainty. We check this expectation when cal- 
culating a lineshape model uncertainty. 



B. Lineshape Model Uncertainties 

We extract f c and v a from their resonance lines by 
binning the excitation attempts into histograms and cal- 
culating the weighted-mean frequency of each line using a 
trapezoid-rule integration. The uncertainty in determin- 
ing f c and v a comes from binomial uncertainties in the 
number of successes in each histogram bin, and we assign 
it as the "statistical" uncertainty. Because the weighted- 
mean values of f c and v a correspond to the same mag- 
netic field — that at the root-mean-square thermal axial 
amplitude in the magnetic bottle — we may use them di- 
rectly in Eq. |10| The assumptions in the weighted-mean 
method are unsaturated lines as well as identical temper- 
ature and drive conditions during cyclotron and anomaly 
excitations. A further assumption is a noise spectrum 
that either is symmetric (and thus does not shift the 
means) or is identical during cyclotron and anomaly ex- 
citations (and thus gives shifts that cancel in the calcu- 
lation of g/2). 

To test these assumptions, we use maximum-likelihood 
fits of the data to a lineshape that includes a specific 
model of the field-noise spectrum: the Brownian-motion 
lineshape convolved with a Gaussian whose width is left 
as a fit parameter. These fits determine the zero-axial- 



amplitude f c and P a , which are used in Eq. 10 to cal- 
culate g/2. The agreement between this line-fit g/2 and 
that from the weighted mean is our primary check on the 
lineshape model, and quantifying this agreement provides 
a systematic "lineshape" uncertainty. 

Specifically, we compare the agreement of two 
weighted-mean and four fitted determinations of g/2. Be- 
cause the weighted-mean method requires binning the 
data into a histogram, we check that our result is in- 
dependent of bin width by using two different numbers 
of bins (50 and 100) at each field. We conduct four 
maximum- likelihood fits to the data: 1. a fit treating 
each excitation attempt separately and fitting the two 
lines sequentially — cyclotron then anomaly; 2. a sequen- 
tial fit to histogrammed data; 3. a simultaneous fit of 
both lines to histogrammed data; and 4. a sequential fit 
to histogrammed data with reduced axial damping. 

In three cases, we use sequential fits of the two lines 
because the cyclotron line better differentiates between 
the two main broadening mechanisms: axial temperature 



and the Gaussian noise function. In the anomaly limit 
of the Brownian-motion lineshape, both broaden the line 
symmetrically; in the cyclotron limit, the Gaussian noise 
function broadens symmetrically while the axial temper- 
ature broadens asymmetrically. We include one fit with 
a reduced axial damping because j z is difficult to resolve 
precisely. Each fit contains six free parameters: the cy- 
clotron and anomaly frequencies and Rabi frequencies, 
the axial temperature, and the width of the convolved 
Gaussian. 

At each field, we treat the 50-bin weighted-mean re- 
sult as the measurement of g/2 and its uncertainty as 
the statistical uncertainty. The maximum one-standard- 
deviation discrepancy, beyond this statistical uncer- 
tainty, of the other determinations is the lineshape un- 
certainty. To be cautious, we avoid a reduction of the 
lineshape uncertainty when averaging the g/2 measure- 
ments at four magnetic fields by treating the minimum 
discrepancy (that at 149.2 GHz) as a correlated uncer- 
tainty. It corresponds to our best understanding of the 
lineshape model. Any additional discrepancy is added 
as an uncorrelated uncertainty. These uncertainties are 
summarized in Table ITT1 

Figure [9] displays the entire dataset, representing 37 
nighttime runs at four magnetic fields. The data, binned 
into histograms (points), fit well to a convolution (solid 
curve) of a Gaussian resolution function (solid inset 
curve) and a Brownian-motion lineshape (dashed curve) , 
as indicated by one-standard-deviation confidence limits 
for distributions of measurements about the fits (gray 
bands). The vertical lines indicate the weighted-mean 
frequencies and their uncertainties from both statistics 
and lineshape model. An additional probe of the broad- 
ening comes from histograms (inset) of the edge-tracking 
data, which is used for drift normalization as described in 
Sec. |III J| Although the precise distribution of this data 
depends on the details of the edge-tracking procedure, 
simulations of our procedure indicate that it should be 
distributed with a width comparable to — within a factor 
of two of — the Gaussian broadening width. As shown in 
Fig. |9j the two agree well. 

The resonance lines at 147.5 GHz and 151.3 GHz ap- 
pear much broader than those at the other fields. This 
additional width appears to be stable throughout the 
data runs at each field (otherwise the narrow lines would 



TABLE II. Summary of the lineshape model analysis. All 
uncertainties are in ppt. 

f c / GHz 147.5 149.2 150.3 151.3 



g- value range 0.73 

statistical uncertainty 0.39 

correlated lineshape 

model uncertainty 0.24 
uncorrelated lineshape 

model uncertainty 0.56 



0.29 0.33 
0.17 0.17 

I 



0.24 




0.24 
0.15 



0.45 
0.24 



0.24 
0.30 
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FIG. 9. Quantum-jump spectroscopy lineshapes for cy- 
clotron (left) and anomaly (right) transitions with maximum- 
likelihood fits to broadened lineshape models (solid) and in- 
set resolution functions (solid) and edge-tracking data (his- 
togram). Vertical lines show the 1-a uncertainties for ex- 
tracted resonance frequencies. Corresponding unbroadened 
lineshapes are dashed. Gray bands indicate 1-a confidence 
limits for distributions about broadened fits. All plots share 
the same relative frequency scale. 



be much noisier at high frequencies) but varies between 
fields. It is consistent with a higher axial temperature, as 
indicated by the asymmetric broadening of the cyclotron 
line with a wider exponential tail, but we have not been 
able to identify the procedural differences at these fields. 

These differences motivate our assignment of separate 
lineshape model uncertainties at each field. The agree- 
ment of the weighted-mean and fitted g/2 determinations 
is much better for the two narrower lines than for the 
wider ones, which is not surprising because the wider 
lines rely more on the lineshape model for line-splitting. 

The weighted-mean and line-fit methods should yield 
a g/2 that is independent of axial temperature. Nev- 
ertheless, we check for any systematic trends related to 
axial temperature by taking an additional set of data at 
149.2 GHz with the refrigerator operating at 500 mK 
instead of 100 mK. Whereas the 100 mK data fit to 
T z = 0.23(3) K, the 500 mK data fit to 0.55(2) K, 
consistent with our deliberate heating. The higher- 
temperature weighted-mean calculation has a statistical 
uncertainty of 0.30 ppt, and the maximum-likelihood- 
fit checks give an uncorrelated lineshape model uncer- 
tainty of 0.46 ppt, both larger than those of the lower- 
temperature data at 149.2 GHz in agreement with the 
temperature-uncertainty correlation noted above. In- 
cluding the statistical and uncorrelated lineshape model 



uncertainties, the difference between the 149.2 GHz, 
500 mK g/2 and the 100 mK g/2 is 0.5(6) ppt, which 
is consistent with zero. 



C. Possible Broadening Sources 

What could cause the additional line-broadening? We 
have modeled this effect above as fluctuations in the 
magnetic field. Attributing the line broadening to field 
noise assumes that the fluctuation timescalc is not so 
fast that the noise averages away during an excita- 
tion attempt. The relevant comparison timescale is the 
inverse-linewidth coherence time (200 /xs for the cyclotron 
line and 200 ms for the anomaly line), and any line- 
broadening noise must fluctuate near to or slower than 
these timescales. Noise-broadening from slow fluctua- 
tions is analogous to the exponential limit (7 Z /Aw <C 1) 
of the noise- free cyclotron line, which takes its shape from 
the long axial fluctuation time and the distribution of 
axial energies. Our edge-tracking data provide an upper- 
timescale of minutes for the noise timescale because we 
see no correlation between adjacent edge-tracking points, 
which come at intervals of several minutes. This range 
of allowed timescales constrains the possible fluctuation 
mechanisms. 

We have considered and ruled out several possi- 
ble sources. Section IIIIDI includes several known 
producers of field noise — the local subway, cryostat 
pressure changes, ambient temperature changes, and 
temperature-dependent paramagnetism — and our efforts 
to reduce them. Phase noise on the cyclotron and 
anomaly drives would mimic field fluctuations, though 
the similar broadening observed on each line is not what 
one would expect given the vastly different drive frequen- 
cies. Estimates based on the microwave signal genera- 
tor's specified phase noise and additional noise from an 
ideal multiplier suggest any frequency deviations should 
be over two orders of magnitude below the level required 
to explain the cyclotron broadening. 

Remaining candidates include relaxation of stresses in 
the solenoid windings or vibration of the trap electrodes 
in an inhomogeneous magnetic field. Even with the mag- 
netic field tuned to its specified homogeneity, a 100 /im 
motion of the trap electrodes would cause a 0.1 ppb field 
variation. Such a motion could be caused by vibrations 
that drive the dilution refrigerator like a 2.2-m-long pen- 
dulum. 

One additional source of broadening is a distribution 
of magnetron radii in the —p 2 z/2 portion of the magnetic 
bottle (Eq. 19). Ideally, our sideband cooling procedure 



will produce a thermal distribution of magnetron ener- 
gies with a temperature related to the axial temperature 
by [28] 



(37) 



(The negative sign indicates that the magnetron de- 
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gree of freedom is unbound.) Such an ideal distribution 
would broaden the cyclotron line by v m /vz times the ax- 
ial broadening Aw, which is 10 2 -10 3 times too small to 
account for the observed width. Nevertheless, we have 
checked that the broadening is not due to a thermal dis- 
tribution of magnetron energies by deliberately sideband 
cooling with the amplifiers on (T z = 5 K) before damping 
the axial energy to 230 mK; we found no change in the 
broadening. Sideband cooling with the coherent distri- 
bution of axial states found during self-excitation would 
produce a coherent distribution of magnetron states with 
a width that could explain the broadening, but the elec- 
tron is not self-excited during sideband cooling. Interest- 
ingly, prior studies of sideband cooling could only achieve 
cooling to energies 400 times the ideal limit [25j. If this 
non-ideal energy limit corresponds to a distribution of en- 
ergies as well, it could explain the observed broadening; 
though, it is unclear what mechanism limits the cooling 
and whether it applies to our setup as well. 

Our knowledge of the lineshape model provides the 
largest uncertainty in recent g/2 measurements [TJ [S]. 
Future efforts will focus on explaining or eliminating the 
lineshape broadening. 



V. CAVITY CONTROL 
A. Overview 

The averaging time needed to observed one-quantum 
cyclotron transitions is obtained by inhibiting [1TJ [35] 
the spontaneous emission of synchrotron radiation. Inhi- 
bition by factors of 200 and more is accomplished with 
a cylindrical Penning trap cavity that was invented for 
this purpose [BJ. The cylindrical trap shapes the radi- 
ation field within the interior of a conducting right cir- 
cular cylinder, modifying the density of radiation states 
in a way that can be studied and understood, at the 
same time as it provides the high quality electrostatic 
quadrupole potential needed to detect one trapped elec- 
tron with good signal-to- noise [7]. 

The extremely useful inhibition of spontaneous emis- 
sion comes at a cost insofar as coupled oscillators - the 
cyclotron oscillator and the radiation mode oscillators 

- pull each other's frequency. The challenge that thus 
arises is that the interaction of the electron and the cavity 
modes also shifts the cyclotron frequency [U [52] . Typi- 
cal cavity shifts are at the ppt-level in the cyclotron fre- 
quency and ppb-level in the anomaly frequency - large 
enough to unacceptably shift the value of g/2. The cylin- 
drical trap geometry was selected to provide familiar and 
well understood boundary conditions within which the 
properties of the radiation field, and hence the shifts that 
these fields produce, could be calculated and understood. 

The Penning trap establishes the boundary conditions 
of a right circular cylinder for microwaves within the trap 

- determining the electric and magnetic field of the trans- 
verse electric and magnetic modes. Of course, the bound- 



aries are not perfect because the electrodes are deliber- 
ately slit so that sections of the cavity can be separately 
biased trap electrodes, because the electrodes contract as 
the trap is cooled from 300 to 0.1 K, and the electrodes 
are not perfectly machined and aligned. The result is 
that the resonant frequency and damping factor for each 
radiation mode are slightly shifted and must be measured 
within the cold trap cavity. 

We used two independent methods to investigate the 
radiation modes of the cylindrical trap cavity in situ. For 
our 2006 measurement, the synchronization of the collec- 
tive motion of many electrons [HI IS] was used to trace out 
the cyclotron damping rate as a function of the electron 
cyclotron frequency, allowing us to identify and label the 
transverse electric and transverse magnetic modes. For 
the 2008 measurement, the one-electron cyclotron damp- 
ing rate was directly measured as a function of both the 
cyclotron frequency and the amplitude of the axial os- 
cillation through the standing wave field of the cavity 
modes. We demonstrate a 165(4) ^m axial offset be- 
tween the electrostatic center of the trap and the center 
as defined by the standing wave fields, and place a limit 
of p < 10 //m on any radial offset. 

The electric and magnetic fields of these modes are 
used as input for an analytic calculation that is prop- 
erly renormalized to avoid self-energy infinities [10L I53j . 
The calculation is then adapted in a semi-empirical way 
to better describe the way that the electron cyclotron 
frequency is shifted and damped. For the 2006 measure- 
ment we were able to reduce the cavity shift uncertainty 
to 0.39 ppt pp. For the 2008 measurement the cavity 
shift uncertainty was reduced by an additional factor of 
six, to 0.06 ppt 0. 



B. Mode detection with synchronized electrons 

Our first technique for probing the radiation modes 
of the not-quite-ideal trap cavity uses the synchronized 
axial motion of approximately 2 x 10 4 electrons trapped 
near its center. When we first developed this method [5J , 
we were able to use it to identify more than 100 radia- 
tion mode of a cylindrical trap cavity, with frequencies 
between 20 GHz and 160 GHz, with those of an ideal 
right circular cylinder. The cavity used then was nearly 
identical to the one used for the 2006 and 2008 measure- 
ments. The method seemed ideal in that it produced 
in situ a spectrum with peaks at the frequency of ra- 



diation modes (Fig. 10), and with lineshapes that were 



Lorentzian and independent of the number of electrons 
used (as long as saturation and strongly coupled regimes 
were avoided). We thus interpreted these widths as the 
inverse of the mode quality or Q factors. 

The electrons are excited by modulating the potential 
applied to the bottom endcap electrode at a frequency 
that is nearly twice the axial oscillation frequency of the 
electrons, 2v z . This parametric modulation of the trap- 
ping potential heats the electrons. The axial oscillation of 
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their center-of-mass grows exponentially but this growth 
is limited to a steady-state value because the trapping 
potential is anharmonic for large oscillation amplitudes. 
Our prior work [51 [9] showed that the measured oscil- 
lation amplitude was large when the cyclotron damping 
rate was high - when the cyclotron frequency of the elec- 
trons was resonant with a radiation mode in Fig. |10| The 
shape of the spectrum was remarkably independent of the 
number of electrons N, which could be varied to change 
the cyclotron damping rate of the center-of-mass motion 
Nj c . In the strong coupling regime, where Nj c exceeded 
the damping rate of the cavity modes, the observed peaks 
split into two. The axial motion of the electrons also gen- 
erates motional sidebands in the observed cavity spectra. 
The method was very robust and reproducible. 

We learned a great deal about how this method worked 
by varying the number of electrons used from 10 5 elec- 
trons [SI H] down to only 2 [S3] , and by varying electro- 
static anharmonicity, the parametric drive strength and 
the slow rate at which the cyclotron frequency swept 
through resonance with the radiation modes. For elec- 
trons driven at twice their preferred oscillation frequency, 
the resulting axial oscillations can have either of two os- 
cillation phases that differ by tt. Between resonances with 
the cavity radiation modes, where the cyclotron damp- 
ing was very small, the electrons oscillated in equal num- 
bers with both of the two axial oscillation phases, so that 
the detected signal from the axial center-of-mass is very 
small. On resonance with the cavity radiation modes, 
with very strong cyclotron center-of-mass damping, the 
electrons synchronize into axial oscillations with predom- 
inately one of the two possible phases, producing a large 
detected signal. 

This method is not yet fully understood despite bound- 
ary conditions that are very carefully controlled. We have 
written down what we believe are the full equations of 
motion for the iV electrons [S], but have not obtained 
a solution that represents the simple and striking behav- 
ior that is observed. One additional clue may be that the 
method worked robustly over a long time in an apparatus 
for which the trap's vacuum enclosure and the detection 
electronics were submerged in liquid helium. In a brief 
trial, we did not succeed in getting robust performance in 
an apparatus in which the trap vacuum container and the 
detection electronics were cooled by thermal contact to a 
liquid helium dewar. Perhaps the temperature and noise 
of the FET detectors are more important than initially 
supposed. 

In the apparatus used for the g/2 measurements, 
cooled to 0.1 K rather than to 4.2 K and with improved 
detection electronics, the interpretation of these cavity 
mode spectra has become much less clear. For causes 
that we have not yet discerned, we sometimes measure 
spectra that are like what we measured in the first studies 



of the cavity spectra [Hj (blue traces in Fig. 10 ). However, 
much more often we measure severely broadened spectra 



(red traces in Fig. 10 1. 

The observed cavity spectra nonetheless give us an im- 
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FIG. 10. Modes of the trap cavity as observed with syn- 
chronized electrons in three separate sweeps. The calibration 
frequencies are indicated with symbols above the plots. The 
striking difference between Maps B and C and Map A is dis- 
cussed in the text. 



portant in situ probe into the spectrum of the cavity ra- 
diation modes. It takes several hours to measure a spec- 
trum. In light of the striking difference (mentioned above 
and discussed in more detail below) the broadened spec- 
tra marked A in Fig. [10] is not used in the cavity shift 
analysis. The 2008 measurements uses the same trap 
cavity and measured spectra as in 2006 [T], but our un- 
derstanding has improved due to the addition of a new 
probe of the cavity spectra (next section) and through 
more thorough analysis. 

The frequency axis for the cavity spectra is calibrated 
during the sweep by applying a strong cyclotron drive 
at discrete frequencies (symbols in Fig. ] 
nant with the magnetic-bottle coupling shifts v 



When reso- 
and 



the parametric drive is no longer resonant, producing 
a sharp dip in the signal such as that on TE241 near 
139.6 GHz. The cyclotron frequencies for the remainder 
of the map are assigned with a linear interpolation func- 
tion between these discrete calibration points. (Map A is 
calibrated with an earlier technique that involved stop- 
ping the magnetic field where marked and measuring D c 
directly.) Comparing the mode frequencies with those of 
an ideal cylinder with dimensions similar to the trap di- 
mensions allows us to identify and label the TE and TM 
modes. 

Three features aid the interpretation of the cavity spec- 
tra 0. 

1. Strong cloud-mode coupling can split the 
Lorentzian response into a pair of normal modes. 
The Q of the split peaks is higher than that of the 
mode, making it difficult to estimate the mode Q 
from parametric mode maps alone. 

2. The large axial motion of the cloud during the 
measurement (limited only by trap anharmonicity) 
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probes the electric-field profiles of the modes and 
amplitude- modulates j c , producing axial-frequency 
sidebands at ±v z for modes with a node at the trap 
center (e.g. TE144 at 145.0 GHz) and ±2D Z for 
modes with an antinode at the trap center. These 
motional sidebands are further explored below for 
the single-electron case. 

3. Non-zero cloud size or a relative offset between the 
electrostatic center and the mode center allows cou- 
pling to modes with nodes at the mode center (e.g. 
TEi 36 at 146.4 GHz and TE 243 at 149.7 GHz). 

For the four magnetic fields at which we measure g/2, 
the modes with the largest influence on the cavity shift 
are TE127, TEi 36 , and TM143. Table |Hl| summarizes the 
mode frequencies and Q values obtained from Lorentzian 



fits to the cavity spectra B and C of Fig. 10 after account- 
ing for any normal-mode splitting. The listed uncertain- 
ties are what is estimated as part of the nonlinear least 
squares fitting to Lorentzian peaks, and does not include 
any additional systematics contribution. 



TABLE III. Comparison of the mode parameters from the 
multi-electron parametric mode maps and the single-electron 
lifetime fits. 
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In light of the differences between the cavity spectra 
measured when this method was first developed, and 
what we now are able to observe, and because we are not 
currently always able to robustly produce such spectra, 
we have carefully reflected upon our measurement tech- 
niqe. Early studies on mode-detection with synchronized 
electrons [9] found the CM amplitude closely followed 7 C 
for wide ranges of parameters, yielding Lorentzian mode 
profiles except in cases of normal-mode splitting or mo- 
tional sidebands. With our current apparatus we find 
three classes of behavior: 

1. Convincing cavity spectra with the detected center- 
of-mass signal tracking 7 C with Lorentzian mode 
profiles and no center-of-mass motion far from 
modes, as seen in B and C of Fig. 
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Broadend cavity spectra in which the signal never 
disappears but increases and decreases with 7 C , as 
seen in A. 

No cavity spectrum at all. With nominally the 
same parameters as for the previous two cases we do 
not see the detected signal change as the magnetic 
field is varied to sweep the cyclotron frequency. 
This is the most common behavior. 



There are substantial differences between the earlier 
apparatus (Ref. [9]) and that used for the 2006-2008 mea- 
surement, even though the trap electrodes themselves are 
essentially the same. 

1. Forty-times lower cavity temperature should in- 
crease the inter-particle Coulomb interaction and 
enhance collective motion [55] |5J>] 

2. Ten-times deeper axial potential should have a sim- 
ilar effect. 

3. Improved detection electronics [35] should offer 
greater detection sensitivity while producing less 
noise to heat the trapped electrons. 

4. Ten-times larger electron clouds are used because 
we find that smaller clouds rapidly saturate, with 
all electrons in one of the two oscillation phases. 

5. More heavily filtered and noise- free electrical en- 
vironment should reduce noise driven transitions 
between the bistable states. 

The three maps of Fig. [10] were taken using the same 
trap cavity over the course of 18 months, during which 
time the electrodes were thermally cycled to room tem- 
perature several times, the refrigerator was inserted and 
removed, and the magnet was quenched with the elec- 
trodes inside. At no point were the electrodes them- 
selves disassembled or adjusted. The general alignment 
of the features and the precise alignment of the calibra- 
tion points in Fig. [10] indicate that the trap cavity and its 
resonant modes are robust against stresses and thermal 
cycles. Any misalignments in the location of a particu- 
lar mode may be attributed to the calibration process, 
specifically to nonlinear charging rates from the power 
supplies, rather than to real shifts in the mode frequen- 
cies. This consistency suggests that the cavity itself is 
stable, though the variety of parametric behaviors dis- 
cussed in the previous paragraph motivates our use of an 
independent, one-electron mode detection technique. 



C. Cavity Coupling to a Single Electron 

1. Overview 

The expressions describing the coupling of a single elec- 
tron to the electromagnetic modes lie at the heart of our 
cavity analysis. Their importance is twofold: first, they 
are required for the calculation of the cyclotron frequency 
shifts, Aw c , and second, measurements of the cyclotron 
damping rate, 7 C , allow a characterization of the cavity 
mode structure independent of the multi-electron tech- 
nique of the previous section. In this section, we present 
two formulas for the cavity-induced cyclotron frequency 
shift and damping rate. They differ in their treatment 
of the electron's axial motion. Equation [40] applies for 
an electron with negligible axial motion, corresponding 
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to the cold thermal axial distribution during cyclotron 
or anomaly excitation. Equation 58 applies for a larger 



axial motion that begins to probe the mode's standing 
wave, corresponding to the self-excited axial state during 
detection and measurement of j c . 



2. Single-mode approximation 

Before beginning a full calculation of the cyclotron 
motion-cavity coupling, it is worth modeling the interac- 
tion between the electron and a single nearby mode, here 
denoted M, to give an indication of the character of the 
electron-mode coupling. This approximation will eventu- 
ally aid in modeling the coupling including axial motion 
in Sec. |V C 5| The interaction may be approximated as 
that of two coupled oscillators with the resulting electron 
frequency shift and damping rate given by |10j 



Auj c 



7c 



7m 



2 1 + S 2 
1 

7M: 



1 + S 2 ' 



(38a) 
(38b) 



Here, 7m is the cyclotron damping rate when the electron 
is exactly resonant with the mode and 5 is the relative 
detuning, defined as 



com 



r M /2 



(39) 



The mode full- width at half-maximum, Tm, arises be- 
cause of losses in the cavity and may be written in terms 
of a quality factor, Qm, with the definition: Qm = 
wm/Tm- The cyclotron frequency is maximally-shifted 
by ±7m/4 at S = ±1. Furthermore, provided the cy- 
clotron frequency is detuned far enough from a mode 
that S > 1, i.e., (w c - w M )/w M > V(2Qm), the shift 
Aw c is Q-independent. 

It is tempting to expand on the single-mode approx- 
imation above by adding the contributions of many 
modes. This mode-sum approach is fundamentally flawed 
because the real part is infinite |53j . A linear divergence 
arises from the inclusion of the electron self-field contri- 
bution to the cavity radiation rather than only the field 
reflected from the walls. A calculation that explicitly 
removes the electron self-field from the cavity standing 
wave, i.e., "renormalizes" the field, yields a finite result 
and is the subject of the next section. 



3. Renormalized calculation of cyclotron-cavity coupling 

While it is possible to tackle the full cylindrical cav- 
ity directly, removing the electron self-field from such 
a calculation is difficult. It simplifies when first ana- 
lyzing the interaction of the electron with two parallel 
plates then adding the contribution of interactions with 
the cylindrical wall. The calculation is quite involved. It 



is presented for a centered particle (z, p — 0) in Ref. [53] 
We extend the calculation to any position in the cav- 
ity. This extension is required for two reasons: our elec- 
trostatic quadrupole suspends the electron slightly off- 
set axially from the cavity mode center (see Sec. VD I, 
and the electron's axial motion modulates the electron- 
cavity coupling when measuring the cyclotron damping 
rate (see Sec. VC5). The derivation for arbitrary posi- 
tion is nearly identical to that in Ref. 53J except we keep 
the terms that vanish for z,p = 0. Rather than rehash 
the lengthy calculation, we state the results and discuss 
some important characteristics. 

For two parallel conducting plates, i.e., the cylindrical 
cavity with po — > oo, the boundary conditions may be 
satisfied with a series of image charges. Renormalizing 
this sum is trivial — simply omit the electron's contribu- 
tion and leave that of the image charges, a result we de- 
scribe below and call Sp. Proceeding to the calculation 
of the full cylindrical cavity and omitting the contribu- 
tion from the endcaps leaves only the correction from the 
cylindrical wall, a result we call Eg. The final result is 
the sum of the contributions from the endcaps and the 
wall. 

At a cavity-shifted cyclotron frequency uj, the fre- 
quency shift Alu c — uj — Q c and damping rate j c is given 
in terms of these two contributions (£_p, £5), the free 
space damping rate 7 c o, and a quality factor Q for all 
modes: 



7c = 



{Ss[( 



;7co 



2Q 



(40) 
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(1 



2Q 



)uj, z 



The correction to g/2 (Eq. 10) is equal to the relative 
shift of the cyclotron frequency: 



Ag cav Aw c 



(41) 



We use this formula in a slightly modified form to calcu- 
late the cavity shifts of the cyclotron frequency (see the 
discussion in the next subsection). 

The renormalized calculation begins by modeling the 
effect of the cavity on the electron as an electric field 
E'(r) arising from image charges in the walls. It modifies 
the transverse equation of motion to read 



— W(r) , 
m '2 



7c0 v = -E'(r). 



(42) 



The longitudinal part of E'(r) gives a negligible correc- 
tion to the trapping potential V(r) [53], but the trans- 
verse part generates the anticipated effects. Using the ra- 
diation gauge, V ■ A = 0, the electric field may be written 
as the time derivative of the vector potential. This vec- 
tor potential satisfies the wave equation with a transverse 
current source and thus may be written as the convolu- 
tion of that source — the moving electron — and a Green's 
function subject to the appropriate boundary conditions, 
see e.g., [33] Sec. 6.3-6.4]. Combining the two transverse 
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velocity components as v 



IVy 



v e~ 



then write Eq. 42 as Eq. 40 where Ep and Eg are pro- 
portional to Fourier transforms of the part of the Green's 
function that arises due to the presence of the cavity 
walls. Note that Ep and Eg are in general complex, 



with the real portion corresponding to a frequency shift 
and the imaginary portion to a modified damping rate. 

As mentioned above, the method of images gives the 
parallel-plate contribution to the cyclotron frequency 
shift and damping rate. As a function of axial position z 
and cyclotron frequency u), it is 



Ep(w, z) = -r 



(43) 



where ro = e 2 /(47reoTOc 2 ) is the classical electron radius 
and F(z) is the Fourier transform of the aforementioned 
Green's function a distance \z\ from an electron or image 
charge: 



F[z) 



1 

R 
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The j = term has been removed from the first sum; this 
exclusion of the electron self-field is the explicit renormal- 
ization required to avoid an infinite result. 

E p depends on an axial offset but not a radial one be- 



cause of the transverse symmetry of two parallel plates. 
Its imaginary part — cyclotron damping — has a sawtooth 
form with sharp teeth where the frequency corresponds 
to an integral number of half-wavelengths between the 
two endcaps (only the odd integers for z = since 
the even integers have a node there). The real part — 
cyclotron frequency shifts — shows peaks at similar inter- 
vals. The nearest such frequency corresponds to eight 
half-wavelengths at 154.5 GHz, far enough away that the 
parallel-plate contribution is smooth in our experimental 
region of interest. 

The contribution from the cylindrical wall is 
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(46) 
(47) 

(48) 



The sums include modified Bessel functions of the first 
and second kinds, 

I v {x)=i- v J v {ix) (49) 
K v (x) = (tt/2) [(/_„(!) - J„(z))/sm(iV)] (50) 

as well as their derivatives. The K' m (p,pPo)/I' m (p p po) 
term comes from the boundary conditions of the TE 
modes, while the K m (p p po) / I m (PpPo) term comes from 



r 



the TM modes. For p — > 0, the Ri functions all go to 
zero except when to = 1, when it goes to 1/2. For z — > 0, 
only the odd-p terms survive. Combined, these limits 
reproduce Eq. 4.28 of Ref. E2 

The frequency shifts and damping from individual 
modes can be seen by looking at the Bessel functions 
in the denominators. For a given p, an increasing w will 
eventually cross a threshold at which p, p becomes zero 
and then imaginary. At that point, we may use the defi- 
nition of I m (x) to substitute 



I m {PpPo)=i m Jrn(P>pPo), 



(51) 



where p, p is the now-real quantity ip p . Since J m (x) and 
J' m (x) have a number of zeros, after uj exceeds the p th 
threshold the sum has poles that may be approximated 



E s (w,z,p) 



\2 

^mnp 



(52) 
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For TE modes, the poles occur when J^(/x p p ) has a 
zero; for TM modes, when J' m (ppPo) has a zero. That is, 
there are poles when u) = io mnp of Eq. 13 Expanding the 
Bessel functions about their zeros yields mode coupling 
strengths 

m "P ~~^2" T,,( x , )J ( x f ) (° 6a ) 

xsin 2 ( ^ ( ^ + 1))jRj(m;a; ; m iL ) 



2rgc 2 1 + sgn(m) / pir 



TM- A 2 



(53b) 



where Xmn are the previously mentioned zeros of the 
Bessel functions and their derivatives. The entire ra- 
dial dependence of the coupling is contained in the Rj 
function, defined by 



Rj(m;x) = -^rJ m (xf - J' m {x) 2 



For zero radius, Rj equals 1/2 if m — 1 and zero other- 
wise. 

Given the above, we can see that the summation in- 



dices p and m in Eq. 45 correspond directly to those in 
the mode indices mnp, and the addition of the m,p th 
term of the sums adds the contributions from all modes 
of that m and p. The threshold above which \i p is imagi- 
nary corresponds to the frequency whose half-wavelength 
fits between the endcaps p times. The single-mode ap- 
proximation of Eq. [52] shows the interaction between the 
electron and a single mode is that of two weakly coupled 
oscillators. 



4- Using the renormalized calculation 



The combination of Ep and Ss in Eq. 40 is the result 
of the renormalized calculation. There, we have included 
cavity dissipation in the from of a mode Q with the re- 
placement io (1 + i/(2Q))io. It is possible to include 
different quality factors for the TE and TM mode classes 
by using Qte in the denominator functions I' m {[i p po) and 
Qtm everywhere else [S3]. It is not possible to include a 
Q for each mode separately. 

The strength of the renormalized calculation is its re- 
moval of the electron self-energy. It has an important 
drawback in that the entire calculation has only four in- 
put parameters: po, zq, Qtm, and Qte- It does not allow 
the input of arbitrary mode frequencies and Qs. If the 
dominant mode-couplings are to one TE and one TM 
mode, then the two mode frequencies and Qs can deter- 
mine the four input parameters. The addition of a third 
mode, however, over-constrains the problem; unless the 
three modes happen to have frequencies that correspond 
to those of an ideal cavity and two happen to share Qs, 
the renormalized calculation will give an incorrect result. 



In this measurement, three modes have the largest in- 
fluence on the electron-cavity coupling: TE127, TM143, 
and TE 136 . Because we can identify the terms in Eq. |45| 
that correspond to a given mode, we modify the renor- 
malized calculation to better approximate our observed 
mode structure. Since the electron is close to (but not 
precisely in) the mode center, the two modes with antin- 
odes at the center (odd p) dominate the coupling, and 
we set the four input parameters of the renormalized 
calculation with the frequencies and Qs of TE127 and 
TM143. These input parameters give an ideal frequency 
and Q for TE136, which we correct to the observed val- 
ues by subtracting the term in Eq. [45] that includes 
the contribution from TE136, — {2tq/zq) sin 2 (37r(z/zo + 
1 ))[-^i(M6Po)/^i(M6Po)] J R/(l;M6P), and adding it back 
with the observed values. Although this moves all 
modes with TEi„6 the next-nearest ones are far away: 
wtei26/(2tt) w 129 GHz and wteW^tt) Rj 169 GHz, 
while wtei36/(27t) Rj 146 GHz. If this were a concern, 
we could do a similar subtraction with the single-mode 

but at the cost of artifacts of 



(54) approximation of Eq. 52 



uncanceled higher-order terms. 



5. Coupling modulated by axial motion 

When the SEO is running, the axial motion through 
the cavity mode field will modulate the coupling at 10 z - 
This motion is present, for example, during cyclotron or 
spin state-detection; it is not present during cyclotron 
or anomaly excitation, the relevant period for any sys- 
tematic shift in g/2. Since measurements of the cy- 
clotron damping rate involve continuous detection of the 
cyclotron state waiting for a decay, we must account for 
this modulation and any amplitude-dependence of the 
damping rate when using measurements of 7 C to deter- 
mine the input parameters to the g/2 cavity corrections. 
This accounting is a nontrivial task, which is intractable 
for the full renormalized calculation but doable for the 
single-mode coupling of Eq. [52] provided that the axial 
amplitude A is much lower than a quarter-wavelength of 
the mode's axial standing wave (A -C zo/p), a condition 
met for typical A. 

All of the z-dependence in the mode-coupling param- 
eters A„ 



x mnp 



of Eq. 53 comes in a single sine function: 



A 



mnp 



sm(f ( 



1))A 



M ■ 



(55) 



where we define Am to be the non-z-dependent part of 
the coupling. For an electron offset z from the center of 
the modes and oscillating at frequency io z with amplitude 
A <C zq/p, we may expand the mode axial-dependence in 
terms of axial harmonics, 

00 

^ f(£± ^M + 1)} A) cosCMt), (56) 

where the fj(z, A) are functions of the axial offset and 
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FIG. 11. Measurement of the cyclotron damping rate at 
146.70 GHz, near the upper sideband of TE136. The cyclotron 
damping rate as a function of axial amplitude (c) extrapolates 
to the desired lifetime. Each point in (c) consists of a damp- 
ing rate measured from a fit to a histogram of cyclotron jump 
lengths (a) as well as an axial amplitude measured from a 
driven cyclotron line (b). 



amplitude. The hrst three are 
f (z,A) = S m(^(f o +l)) 

/ 1 ( z ,A) = cos(^(^ + l)) 
/ 2 (z,A) = sin(f (£ + 1)) 
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Including this expansion of the axial oscillation in 
the transverse equation of motion [571 App. A] yields 
an amplitude-dependence to the single-mode coupling 
strength as well as a series of axial harmonics to the mode 
frequency, lom- 



Ac 



■ 1c A^w 



1 



(58) 



UJ 2 - (0J M - ju z ) 2 UJ 2 - (U)M + ju z ) 2 



As before, we may include a damping width by substi- 
tuting uj — > (1 + i I (2Q))u) in the two fractions within the 



usual single-mode coupling of Eq. 52 

D. Single-electron mode detection 

Using the above descriptions of electron-cavity cou- 
pling, we measure the cyclotron damping rate as a func- 
tion of cyclotron frequency and of position in the trap to 
determine the location and Q of the three closest coupled 
modes and to characterize the alignment of the electro- 
static and mode centers. Since we are able to perform 
a QND measurement on the cyclotron state, measuring 
the cyclotron damping rate simply consists of making 
many (typically hundreds of) jumps and fitting the dis- 
tribution of jump lengths to a de caying exponential with 
time-constant 7" , as in Fig. 



Hi. Because the axial mo- 



tion is self-excited during detection, we must measure the 
damping rate as a function of amplitude. The amplitude- 
dependence goes as even powers of A (Eqs. 57 & 58 1, 
and since the amplitude is much less than a quarter- 
wavelength for the relevant modes, terms of higher-order 
in A get progressively smaller, allowing approximation as 
a quadratic function, 



j{A) = 70 + 72A 2 



(59) 



We measure the axial amplitude with the driven cy- 
clotron lineshape of Eq. 30 where the driven axial mo- 



tion in the magnetic bottle causes the electron to see a 
higher average magnetic field, resulting in a cyclotron fre- 
quency shift as in Fig. [Tip. Figure [TTj: shows an example 
of the measured damping rate as a function of amplitude 
close to the upper axial sideband of TE136 (^tei36 + v z ). 
It displays a large amplitude-dependence in j c because 
of the proximity of v c and the sideband, which becomes 
more prominent as larger axial oscillations increase the 
modulation of the mode coupling. 

After repeating such measurements at many cyclotron 
frequencies, we amass two sets of data, the zero- 
amplitude cyclotron decay rates 70 and the quadratic 
amplitude-dependence coefficients 72, as functions of D c . 
From these data, we extract the frequencies of the three 
nearest coupled modes: TE127, TEi 36 , and TM143. Since 
the lifetimes are heavily Q-dependent, we must include 
three Qs, bringing the total number of fit parameters to 
six. 

The renormalized calculation provides the model for 
fitting the zero-amplitude cyclotron decay rates. As de- 
scribed above, we use its four free parameters to deter- 
mine the frequencies and Qs of TE127 and TM143 and add 
the TE136 parameters by subtracting the TEi„g term in 
the renormalized calculation for the calculated frequency 
and Q and adding it back in with the frequency and Q 
as fit parameters. 

For the quadratic-amplitude-dependence data, we use 
the A 2 term in the single-mode coupling expansion of 
Eq. [58] to write 
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FIG. 12. Measurement of the offset between the electrostatic 
and mode centers. Here, z refers to the position relative to 
the electrostatic center, and the minimum in the cyclotron 
damping rate corresponds to the mode center. 



where the sum is over the three modes of interest. 

We form a x^ using the zero-amplitude data and the 
renormalized calculation and a x^ 2 using the quadratic- 
amplitude-dependence data and Eq. [60] Since the elec- 
tron is close to centered axially (we discuss the offset be- 
low and include the measured offset in the fits) , the zero- 
amplitude data is more sensitive to the two modes with 
central antinodes (TE127 and TM143) and the amplitude- 
dependence data to the sidebands of TE 136 , which has a 
central node. The fit consists of minimizing the two x 2 s; 
though the proper weighting of the two is not clear a 
priori, it makes little difference to the results. 

Figure [T3pfec display the lifetime data and fits, and Ta- 



ble III lists the results and compares them to those from 
the parametric mode maps. The two independent meth- 
ods should agree but do not. When calculating the cav- 
ity shifts, we assign uncertainties large enough to include 
both results for the mode frequencies (see Table [rv] ) . For 
the mode Qs, to which the cavity shifts are largely insen- 
sitive, we use the results from the lifetime fits because 
of the strong-coupling ambiguity in the parametric mode 
map values. 



Electron position in the cavity 

Knowledge of the electron position relative to the cav- 
ity modes is important for calculating the electron-mode 
coupling and thus the cyclotron frequency shift. We 
determine the axial misalignment by measuring the cy- 
clotron damping rate as a function of z, using antisym- 
metric endcap potentials to move the electron along the 



trap axis [BJ. Figure 12 plots such a measurement with 
the cyclotron frequency tuned midway between TEi3g 
and its upper sideband — close enough to TE136 to give a 
large ^-dependence but far enough detuned that the cy- 
clotron lifetime is long enough to see single excitations. 
Since the mode-couplings are even functions of z, the 
damping rate should go as z 2 with an extremum at z = 
for the modes. Figure fl2] shows this extremum, which fits 
to a misalignment of 165(4) /im. In addition, since the 
extremum is a minimum, the nearest coupled mode must 
have a node at the mode center (even p) , demonstrating 
that the TE127/TE136 identification in Fig. 10 is correct. 



The cause of this offset is not known, though the con- 
sistent presence of TE 136 in the parametric mode maps 
suggests that it is stable. Measurement of a third trap 
"center," the minimum of the magnetic bottle, agrees 
with the electrostatic center. Attempts to model the ob- 
served offset as asymmetric spacing between electrodes 
does not yield a convincing explanation and more exotic 
trap deformations such as a tilted endcap or a compensa- 
tion electrode that protrudes slightly into the cavity are 
difficult to model. Because both the parametric mode 
maps and Fig. 12 indicate an offset, we include it when 

We build our confidence 
at four cyclotron 



calculating the cavity shifts, 
in this procedure by measuring g/2 
frequencies with different cavity shifts and showing the 
agreement between the predicted and measured shifts. 



We estimate the radial alignment of a single electron by 
tuning its cyclotron frequency into resonance with three 
modes that have nodes at the radial center, i.e., have 
m. ^ 1, and comparing the measured cyclotron damping 
rate to that predicted by the renormalized model with 
p — 0. Since the m ^ 1 modes do not couple to a 
radially-centered electron, a measured cyclotron damping 
rate that is faster than the calculated damping rate could 
indicate a radial misalignment. For the cases where we 
observe such a discrepancy, we use the full, p-dependent 
renormalized calculation to estimate the range of radial 
offsets that could explain the observed damping rates 
(all calculations include the axial offset of the previous 
paragraphs). In each of the three cases (TE 03 5, TM027, 
TE 4 3 ), we measure damping rates close to that predicted 
for p — 0, and we set the limit p < 10 /im. 
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E. Cavity-shift results 

We calculate the cavity shifts and uncertainties from 
the mode parameters and their uncertainties (Table IV) 
via the renormalized calculation. The cavity shifts are 
independent of mode Q for cyclotron frequencies with 
relative detunings (cD c — wm)/wm ^> 1/(2Qm)- Because 
all our g/2 measurements meet this requirement easily, 
the cavity shift uncertainty comes only from the errors 
in mode frequencies. Because of the axial offset, mode 
TE136 affects the shifts, and we calculate with this mode 
at its observed frequency using the mode-moving tech- 
nique described in Sec. |V C 4| Because fits to the para- 
metric mode maps and to the single-electron lifetime data 
yield slightly different frequencies for the three nearest 
coupled modes (see Table |III[ ), we assign uncertainties 
large enough to include both. The trap-radius limit only 
has a significant effect near two modes, TE243 and TE043, 
and we again use the mode-moving technique to place 
them at their observed frequencies though it makes little 
difference to our result because none of our four measure- 
ments of g/2 were resonant with either mode. Because it 
appears in the parametric mode map, we include TM027 
in the calculation, although it does not change the result 
noticeably. 

Fig ure [T3p lfee display the results of this analysis, and 
Table |VI| shows the calculated cavity shifts for our four 
measurements of g/2. The shifts span over 10 ppt with 
uncertainties around 100 times smaller than that range. 
The lowest uncertainties are below a part in 10 13 , over 
six times smaller than in our 2006 measurement [T] and 
low enough that this systematic uncertainty is no longer 
a dominant error in the measurement of g/2. 



VI. POWER SHIFTS 
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FIG. 13. Cavity shift results come from synchronized elec- 
trons (a) and from direct measurements with one electron of 
7 C (b) and its dependence on axial amplitude (c). Together, 
they provide uncertainties in the frequencies of coupled cavity 
modes (gray) that translate into an uncertainty band of cavity 
shifts Ag cav /2 (d) whose half-width, i.e., the cavity shift un- 
certainty, is plotted in (e). The diamonds at the top indicate 
the cyclotron frequencies of the four g/2 measurements. 



We expect neither v a nor f c to shift with cyclotron 
or anomaly power, but previous measurements of the 
electron g/2 at the UW showed unexplained systematic 
shifts of the cyclotron frequency with both drive pow- 
ers [H . The origin of these power shifts in the UW 



TABLE IV. Parameters used in calculating the cavity shifts. 
For comparison, we include earlier estimates from the 2006 
measurement, which used the same trap cavity. 







2006 


this 






measurement [1] 


measurement 


TE127 


v c / GHz 

Q 


146.350(200) 
> 500 


146.309(27) 
4900(300) 


TE136 


v c 1 GHz 




146.428(15) 


Q 




4800(200) 


TM 143 


v c 1 GHz 

Q 


151.900(200) 
> 500 


151.832(37) 
1270(70) 


electrostatic 
offset 


z / fim 
p/ fim 






165(4) 
< 10 



measurements remains unknown |29j . and extrapolation 
to zero cyclotron power involved correcting shifts of sev- 
eral ppt in g/2 [4j 129) . Estimates comparing our drive 
power to that used in the UW measurements suggest that 
our narrower lines and single-quantum cyclotron tech- 
nique require drive strengths low enough that the power 
shifts are negligible. We have yet to see a power shift in 
our apparatus, though experimental searches are time- 
consuming and the statistical uncertainty in the current 
search is comparable to our final uncertainty in g/2. 



A. Anomaly power shift estimate 

The UW experiment showed anomaly power shifts 
of several ppb in the anomaly frequency [29] . An 
off-resonant anomaly drive during cyclotron excitation 
shifted the cyclotron line by a similar amount, and the 
two shifts canceled in the frequency-ratio calculation 
of g/2 4J. The origin of these shifts is unknown, al- 
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though experiments with a variable-strength magnetic 
bottle showed that they increase with the magnitude of 
the bottle strength, independent of its sign [53] . 

Direct comparisons between the anomaly power used 
in the UW experiment and that used here are difficult 
because the experiments use different anomaly excitation 
techniques. The UW excitations were primarily driven 
with countcrflowing current loops in split compensation 
electrodes, while we drive the electron axially through the 
zpp gradient of the magnetic bottle. Unlike the current- 
loop excitation technique, our axial-excitation technique 
provides a clear mechanism for an anomaly power shift 
by increasing the average axial amplitude and, therefore, 
the average magnetic field seen by the electron (because 
of the z 2 z part of the magnetic bottle). We estimate this 
shift below and expect it to be both smaller than our 
current precision and canceled in the calculation of g/2 
by a similar cyclotron shift from a detuned anomaly drive 
during cyclotron excitation. 

For a driven axial amplitude, z a , the frequency shift 
from the motion through the magnetic bottle is 



B 2 



(61) 



Were this shift to approach the linewidth, we would be 
forced to use the driven lineshape of Eq. [30] when extract- 
ing frequencies. In order to calculate the expected size 
of the power shift, we must estimate z a ; only amplitudes 
over 800 nm will produce shifts at the 0.1 ppb level in 
frequency. We estimate z a using two methods: the ob- 
served anomaly transition rate and a calibration of the 
drive voltage. The estimates give similar amplitudes, and 
neither predicts anomaly power shifts at our precision. 

From Eq. |35[ we see that the anomaly transition rate 
goes as the product of the Rabi frequency squared fi 2 . 
times the lineshape function x( w )- The anomaly Rabi 
frequency goes as z a |50j . and we can estimate z a from a 
typical peak excitation fraction, P p k, using 



P pk = -Tfl aX {oJ p k) 



— T I — BoZn 1 

2 I 2 2m 1 



(62) 



mh(uj c - uj m ) 



X(^pk)- 



Because the lineshape function is normalized to unity, 
its value on-peak in the Lorentzian lineshape limit — the 
limit corresponding to the anomaly line — is inversely- 
proportional to the linewidth 



, s 2 /2Aw 2 
X(wp k ) = h 7, 



(63) 



For typical experimental parameters, we must drive to 
z a w 100 nm to achieve a 20% excitation fraction. 

Alternately, we can estimate the driven amplitude 
based on the rf voltage on the bottom endcap. An end- 
cap driven with amplitude V a excites the electron to an 



amplitude given by [22 
Za ~ 2z 



V* 
Vr 



(64) 



(Vr is the ring electrode potential, and ci, d, and zq are 
geometric factors equal to approximately 0.78, 3.5 mm, 
and 3.8 mm in our trap.) We calibrate the drive ampli- 
tude using the anomaly-power-induced axial frequency 
shifts discussed in Sec. |III| They indicate 30 dB of atten- 
uation between the anomaly frequency synthesizer and 
the trap electrode. This agrees with the attenuation we 
measure in the drive line during room-temperature cali- 
brations and seems reasonable given the 20 dB cold at- 
tenuator installed at the IK pot and some additional loss 
in the stainless steel semi-rigid coaxial cable. The high- 
est anomaly power used for g/2 data was -16 dBV at the 
synthesizer, which would attenuate to V a — 5 mV at the 
bottom endcap and drive the electron to z a — 250 nm. 

Driven axial amplitudes around 100-250 nm should 
only shift the anomaly and cyclotron frequencies at the 
1-10 ppt-level, which is far too small to affect the line- 
shapes (the error in g/2 would be lower by 1000 if the cy- 
clotron and anomaly shifts were uncorrelated, but should 
be even smaller because the shifts cancel in the frequency 
ratio). 



B. Cyclotron power shift estimate 

Unlike the anomaly power shifts, the cyclotron power 
shifts seen in the UW experiments did not cancel in g/2 
and added a 1.3 ppt uncertainty to their 1987 result [4]. 
The shifts appeared as a resonant effect of unknown ori- 
gin with a resonant drive shifting the cyclotron line sev- 
eral ppb but a detuned cyclotron drive with the same 
power not shifting the anomaly line. Investigations in 
a trap with a variable-strength magnetic bottle showed 
that the shift scaled with B2 in magnitude and sign. 
In [29] , the authors hypothesize that the shift could have 
originated in an excitation of the magnetron motion be- 
cause a typical shift could have been explained by a 10% 
increase of the magnetron radius. 

Our cyclotron excitation technique, injecting mi- 
crowaves into the trap cavity, is similar to that used in the 
UW measurements, so we can compare our technique to 
theirs. Our lower temperature narrows the lines by a fac- 
tor of ten, requiring less power to drive transitions. The 
measured bottle-dependence suggests that our ten-times- 
stronger magnetic bottle could cancel the advantage of 
our narrower lines. The overall shifts should still be 
reduced because our single-quantum-jump spectroscopy 
only needs to excite to the n = 1 state less than 20% 
of the time. At the UW, typical excitations sustained 
the electron at energies corresponding to n > 4 [28] [29] . 
Naively, exciting to an average energy of n = 4 requires 
20-times more power than an average energy of n = 0.2, 
and this power reduction alone would reduce several-ppt 
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FIG. 14. The lower plots show cyclotron (a) and anomaly (b) 
data taken in search of cyclotron and anomaly power shifts. 
The upper plots compare the weighted-means of these lines. 



TABLE V. Summary of power-shift searches 



test 


"shift" / ppb 


/ c with double anomaly power 


-0.18(13) 


/ c with half cyclotron power 


0.13(19) 


v a with double anomaly power 


0.11(35) 


v a with double cyclotron power 


0.38(35) 



shifts in g/2 below our precision. The relativistic shifts 
between cyclotron levels suggest additional power in the 
UW drives because excitations above n = 1 involve driv- 
ing in the exponential tails of the higher states' reso- 
nances. In addition, if the power shift is indeed related to 
driving a magnetron-cyclotron sideband, our ten-times- 
higher magnetron frequency and ten-times-narrower cy- 
clotron lines put the closest magnetron sideband, which 
was fewer than 10 linewidths away at the UW, 100 times 
farther from the cyclotron resonance. 



C. Experimental searches for power shifts 

Although we do not expect any cyclotron or anomaly 
power shifts of v a or f Cl their existence in the UW mea- 
surements makes us proceed with caution and look for 
them anyway. We examine the shifts of each line individ- 
ually to ensure that no systematic effects (even ones that 
cancel in g/2) go unnoticed. We look for a cyclotron fre- 
quency shift by running three cyclotron scans: a control, 
one with double the detuned anomaly power, and one 
with half the cyclotron power (lower to avoid saturation). 
The scans are interleaved in the same way we interleave 
cyclotron and anomaly scans during g/2 measurements, 
alternating single sweeps of each line and including edge- 
tracking to remove long-term drifts (see Sec. Ill I. The 
resulting cyclotron lines are shown in Fig. 



14 1. We 



calculate the cyclotron frequency of each line with the 
weighted-mean method (the offset from f c cancels when 
subtracting for a frequency shift). Frequency differences 
between methods are summarized in Table [V] 

To look for anomaly frequency shifts, we run three 
anomaly scans — a control, one with double the detuned 
cyclotron power, and one with double the anomaly power 



(the control power is low enough that we can double the 
power without saturating) — interleaved and normalized 
via edge-tracking as before. The resulting anomaly lines 
are shown in Fig. |14b , which includes the frequencies 
calculated by the weighted-mean method. Table [V] sum- 
marizes the differences. 

The results in Table [V] as consistent with zero. The 
largest "shift" is that of the anomaly frequency with 
cyclotron power — the only one of the four not seen at 
the UW. The data of Table [V] suggest that any power 
shift will be < 0.35 ppb in frequency, which is con- 
sistent with the limits of our prior studies (detailed in 
Sec. 6.2 of Ref. [30]) and with our expectation of no 
shift at our current precision. The uncertainties are lim- 
ited by our ability to resolve the lines in a timely man- 
ner. (The number of nights spent assembling the data in 
Fig. [14] exceeds half the number used to determine the 
g/2 value.) The anomaly line in particular requires the 
time-consuming discrimination between |0,f) and |1,4-) 
after each anomaly pulse, and any search for a system- 
atic shift in the anomaly frequency multiplies the number 
of times this must occur. Because we estimate that no 
power shift should occur and our experimental searches 
are limited by our ability to resolve the lines, we apply 
neither a correction nor any additional uncertainty from 
power shifts. 



VII. RESULTS AND APPLICATIONS 



A. Most Accurate Determination of the Electron g 

Value 



The result for the electron magnetic moment in Bohr 
magnetons, 

g/2 = 1.001 159 652 180 73 (28) [0.28 ppt], (65) 

comes from the weighted average of the four measure- 
ments with uncorrelated and correlated uncertainties 
combined appropriately. The result has 2.7 and 15 times 
lower uncertainty than the 2006 and 1987 measurements 
and 2 300 times lower uncertainty than has been achieved 
for the heavier /i lepton [5 . Table VI summarizes the 



measured values, shifts, and uncertainties for the four 



TABLE VI. Measurements and shifts with uncertainties, mul- 
tiplied by 10 12 . The cavity-shifted "g/2 raw" and corrected 
"g/2" are offset from our result in Eq. 65 



/ c 147.5 GHz 


149.2 GHz 


150.3 GHz 


151.3 GHz 


g/2 raw -5.24(0.39) 


0.31(0.17) 


2.17(0.17) 


5.70(0.24) 


Cav. shift 4.36(0.13) 

Lineshape 
correlated (0.24) 
uncorrelated (0.56) 


-0.16 (0.06) 

(0.24) 
(0.00) 


-2.25 (0.07) 

(0.24) 
(0.15) 


-6.02 (0.28) 

(0.24) 
(0.30) 


g/2 -0.88(0.73) 


0.15 (0.30) 


-0.08 (0.34) 


-0.32 (0.53) 
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FIG. 15. Four measurements of g/2 without (open) and with 
(filled) cavity-shift corrections. The light gray uncertainty 
band shows the average of the corrected data. The dark gray 
band indicates the expected location of the uncorrected data 
given our result in Eq. [65] and including only the cavity-shift 
uncertainty. 



separate measurements of g/2. The uncertainties are 
lower for measurements with smaller cavity shifts and 
narrower linewidths, as might be expected. We no longer 
quote uncertainties for variations of the power of the 
v a and f c drives; although power-shifts appeared in the 
1987 measurement [UHS]; our narrower lines and single- 
quantum cyclotron technique require much lower drive 
powers and we estimate — and check experimentally — 
that they are no longer important. 



Relation of 2006 and 2008 measurements 

The 2008 measurement (SJ is an independent measure- 
ment that is consistent with the 2006 measurement pQ. 
Essentially the same apparatus is used to make a fresh 
data set. However, the apparatus is now better under- 
stood and both the measurement and analysis procedures 
are significantly improved. For example, the electron is 
used as a relative magnetometer to allow many data sets, 
measured on different days, to be combined. This gives 
lineshape curves with a signal-to-noise good enough to 
compare their shape with theoretical expectations. (Pre- 
viously, each day's measurements were combined to get 
a g value, and the values from different days were aver- 
aged.) 

The biggest reduction in the uncertainty in our second 
measurement comes from a better understanding of cav- 
ity shifts. Two independent probes of the cavity mode 
structure allow identification of nearly all modes and the 
quantification of an offset between the effective center of 
the trap for the radiation modes and for the electrostatic 
quadrupolc potential. By measuring g/2 at four mag- 
netic fields with cavity shifts spanning over 30 times our 
final uncertainty, we precisely test this once-dominant 



uncertainty, and demonstrate that we can assign an un- 
certainty that is much smaller than that estimated for 
our first measurement. 

Retroactively applying the improved understanding 
and the modified analysis developed for the second mea- 
surement to the first would require starting from the raw 
data. Improvements in measurement methods cannot be 
retroactively implemented, of course. We thus believe 
that the 2008 measurement should be regarded as super- 
seding the 2006 measurement rather than trying to av- 
erage the two measurements, which would only insignifi- 
cantly change the value of g/2 in the second digit of the 
2008 uncertainty. The correlations between the possible 
systematic uncertainties that limit the two measurements 
has neither been studied nor reported carefully enough to 
allow an appropriate averaging of the two measurements. 



B. Most Accurate Determination of a 

The new measurements of the electron g/2 determine 
the fine structure constant about 12 times more accu- 
rately than the next-most-precise method (Fig. [2| . The 
relationship between g/2 and a has been summarized in 
detail in [22], with the final value updated in [2]. Here 
we give only the results and a brief summary. 

Within the standard model of particle physics the elec- 
tron g/2 is related to a by 



+ G (2) +Cl (2)" + ft (2)' + ft (2 



"T Q-hadronic T" ^-wcak- 



(66) 



The leading contribution to g/2 is the 1 that is predicted 
for a Dirac point particle. Vacuum fluctuations modify 
the interaction of the electron with the magnetic field, 
increasing the effective magnetic moment of the electron 
by approximately one part per thousand. This addition 
is described by the infinite QED series in powers of a/ir, 
with coefficients C n determined by n-vertex QED calcu- 
lations for the interaction of electrons and photons. The 
first three coefficients (C2, C4, and Cq) have all been cal- 
culated exactly. A substantial numerical calculation has 
determined C%, and a numerical calculation of C10 is un- 
derway. A related series involving the \i and r leptons 
yields a small contribution, a^, T . Much smaller hadronic 
and weak contributions, dhadronic an d a W cak> have been 
calculated accurately enough that they do not add un- 
certainty at the current level of precision. References to 
the most recently calculated values are provided in [22] , 
The fine structure constant is determined from the 



measured g/2 by solving Eq. 66 for a to obtain 



a- 1 = 137.035 999 084(33) (39) [0.24 ppb] [0.28 ppb], 
= 137.035 999 084(51) [0.37 ppb]. (67) 

The first line shows experimental (first) and theoretical 
(second) uncertainties that are nearly the same. The 
uncertainty in a is now limited a bit more by the need 
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for a higher-order QED calculation (underway [21]) than 
by the measurement uncertainty in g/2. 

In more detail, the theory uncertainty contribution to 
a is divided as (12) and (37) for C$ and Cio- It should 
decrease when a calculation underway |24| replaces the 
crude estimate C w = 0.0 (4.6) [H1E2]- The a" 1 of Eq.[67 



will then shift by 2a 3 ir- 4 C w , which is 8.0 C w x 10" 9 .^ 
change A 8 in the calculated C$ = —1.9144(35) would 
add 2a 2 7r- 3 A 8 . 

The independent methods for determining a that come 
closet to our accuracy are the "atom-recoil" measure- 
ments, so called because their uncertainty is limited by 
measurements of recoil velocities in Rb and Cs atoms. 
They rely on many experiments, including the measured 
Rydberg constant [551 155]. the Rb or Cs mass in amu [50] . 
and the electron mass in amu [BTJ [55]. The needed 
/i/M[Rb] comes from a measurement of the recoil of a 
Rb atom in an optical lattice (25j [27l [63] • The needed 
h/M[Cs] comes from an optical measurement of the Cs 
Dl line [2BJ and the "preliminary" recoil shift for a Cs 
atom in an atom interferometer [B3] . Although these de- 
terminations of a have an uncertainty that is currently 
12-22 times larger than ours, improvements are expected 
in experiments that are underway [651 166] . 



C. Most Precise Test of QED 

The most stringent test of QED, to the highest order 
in a/'K, comes from comparing the measured g/2 to the 
value calculated using Eq. 66 using the best available 
value of a that is not determined from the electron g. 
Our latest g, compared to Eq. 66 with a(Rb), gives a 
difference gl [27] 



\Sg/2\ < 8 x 10" 



(68) 



The good agreement testifies to the remarkable success of 
QED. The prototype of modern physics theories is thus 
tested far more stringently than its inventors ever envi- 
sioned [57] . 

The latest g/2 measurement is now accurate enough to 
allow a 10-20 times more stringent test of QED, should 
a comparable-accuracy measurement of a become avail- 
able. We thus strongly emphasize the compelling need 
for greatly improved independent measurements of a. 



D. Limits on Electron Substructure 

The same comparison of the measured g/2 and the 
value calculated from Eq. 66 using the best available in- 
dependent a probes the internal structure of the electron 
[221 |B5]. A composite electron is constrained to have con- 
stituents with a mass m* > m/y/Sg/2 = 180 GeV/c 2 , 



available) then we could set a limit m* > 1 TeV. These 
high energy limits seem somewhat remarkable for an ex- 
periment carried out at 100 mK. However, a search for a 
contact interaction in electron-positron collisions at LEP 
sets a more stringent limit, m* > 10 TeV (55) . 



E. Test of CPT Invariance with Leptons 

Already the most precise test of CPT invariance with a 
lepton system comes from comparing the measured mag- 
netic moment of the positron and the electron [4]. A 
new measurement underway at Harvard aims to improve 
the sensitivity of this test by a factor of 15 or more, by 
applying the demonstrated new electron methods to a 
positron. 



F. Application to Dark Matter 



The comparison of the measured g/2 and the value cal- 



corresponding to an electron radius R < 1 x 10 m. 

If this test was limited only by the experimental uncer- 
tainty in g/2 (i.e. if a much better independent a becomes 



culated from Eq. 66 using the best available independent 
a is also relevant to one model that attempts to explain 
dark matter. The measured g/2 is accurate enough to al- 
low the discovery of, or to rule out, proposed dark-matter 
particles with a mass that is close to the electron mass 
[70) . if and when a more accurate independent measure- 
ment of a becomes available. 



VIII. OUTLOOK 

The new g/2 prepares the way for further tests of the 
standard model, pending the availability of an indepen- 
dent a at the uncertainty reported here. In addition, 
the techniques used to measure the electron g/2 clear 
the way for a series of new measurements, some of which 
have already begun. 

First, measuring the positron g/2 to the same preci- 
sion would improve upon the most stringent lepton CPT 
test [4] and constrain possible violations of Lorentz in- 
variance [71] . Except for the loading mechanism and 
an inverted ring voltage, a positron g/2 measurement 
would proceed identically to the electron measurement 
presented here. 

Second, a direct measurement of the proton-to-electron 
mass ratio would combine the sub-ppb electron cyclotron 
frequency resolution presented here with existing tech- 
niques for 90 ppt resolution of the proton cyclotron 
frequency [72 to compete with the existing 0.4 ppb 
limit [20]. 

Third, recent observations with a single trapped pro- 
ton [75] open the way to proposed direct measurements 
of the proton and antiproton magnetic moments at the 
ppb-scale [TTJ[74]. These would reduce the existing uncer- 
tainties by factors of 10 and 10 6 , respectively, and provide 
an important test of CPT invariance [75]. The challenge 
to such a measurement is QND detection of a spin flip 
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in a magnetic bottle because the smaller magnetic mo- 
ment and larger mass reduce the axial frequency shift by 
over 10 4 , though a larger bottle gradient can compensate 
for some of the reduction. Our goal is to realize with a 
proton and antiproton the great signal-to-noise ratio and 
sensitivity to frequency changes that have been realized 
with the one-electron SEO. 

Fourth, access to the lowest quantum states of a 
trapped electron and the lack of radiative damping of 
any degree of freedom except cyclotron motion have led 
to several quantum information proposals using electrons 
in Penning traps, e.g., [7SHT9"] that perhaps could be re- 
alized in a carefully optimized planar Penning trap [SU] . 

IX. CONCLUSION 

In conclusion, precise control of the location of and 
the coupling to the electromagnetic modes of the elec- 
trode cavity reduces the once-dominant cavity shift un- 
certainty. This results in a measurement of the electron 



magnetic moment 15 times more accurate than the 1987 
measurement that provided the best g/2 and a for nearly 
20 years. With the measurement limited by the resolu- 
tion and model of the cyclotron and anomaly lines, fu- 
ture work on the electron g/2 should focus on enhancing 
magnetic field stability, narrowing the lines, and building 
signal-to-noise. The techniques used in this result may 
be directly applied in measuring the positron g/2 and 
may be adapted to a direct proton-to-electron mass ratio, 
to measurements of the proton and antiproton magnetic 
moments, and to trapped-electron quantum information 
studies. With an independent a of similar precision, the 
new g/2 would make possible 10-times more stringent 
tests of extensions to standard model. 
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